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ABSTRACT 

This  study  is  an  attempt  to  reduce  the  effect  of  the  leakage  vortex  in  axial  flow 
turbines.  A  3D  Navier-Stokes  CFD  solver  with  k-s  turbulence  modeling  was  used 
compute  the  flow  through  an  axial  flow  turbine  with  modified  blade  tip  designs.  The 
three  modified  blade  tip  designs  each  involved  adding  a  chamfer  to  the  tip  of  the  blade,  in 
an  attempt  to  diffuse  the  leakage  flow  through  the  gap  and  obstruct  the  leakage  flow  with 
the  outer  casing’s  shear  layer.  A  baseline  flat  tip  case  and  three  modified  tip  cases  were 
simulated  and  the  leakage  flow  and  vortex  for  each  was  analyzed  in  detail.  By 
chamfering  the  blade  tip  near  the  trailing  edge  of  the  gap,  the  leakage  flow  inside  the  gap 
was  turned  toward  the  direction  of  the  blade’s  camber.  This  turning  resulted  in  a 
decrease  in  the  size  and  strength  of  the  leakage  vortex  and  its  subsequent  losses,  while  at 
the  same  time,  did  not  reduce  the  blade  loading  by  an  appreciable  amount. 
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CHAPTER  1:  INTRODUCTION 


The  tip  leakage  flow  is  one  of  the  most  prevalent  and  influential  features  of  the 
flow  through  turbomachine  rotors.  Along  widi  its  resulting  vortex,  the  leakage  flow  acts 
as  a  source  of  both  rmsteadiness  and  aerothermal  loss,  and  is  therefore  detrimental  to 
turbomachine  performance.  Furthermore,  the  leakage  flow  is  a  phenomenon  that  is 
difficult  to  either  predict  or  measure  in  most  turbomachines.  The  research  presented  in 
this  report  is  a  computational  attempt  to  reduce  the  leakage  flow  and  vortex  in  axial  flow 
turbines  through  blade  tip  modification. 

1.1;  A  Description  of  the  Tip  Leakage  Flow  and  its  Environment 

1.1.1:  Physical  Explanation  of  the  Tip  Leakage  Flow 

The  tip  leakage  vortex  is  a  direct  result  of  the  loading  on  the  blade,  and  can  be 
described  two  different  ways: 

Description  1:  At  the  blade  tip,  there  is  an  open  path  for  fluid  to  flow  from  the  pressure 
surface  to  the  suction  surface,  as  well  as  a  pressure  gradient  to  drive  the  flow.  This  is 
known  as  the  tip  leakage  flow,  and  is  illustrated  in  figure  1.1.  The  leakage  jet  interacts 
with  the  main  passage  flow  and  rolls  up  into  a  tip  leakage  vortex.  The  intense  mixing  of 
the  two  streams  (leakage  flow  and  main  passage  flow)  causes  major  viscous  losses. 
Description  2:  As  with  any  cambered  airfoil,  the  velocity  difference  between  the  pressure 
and  suction  surfaces  results  in  a  bounded  circulation  around  the  blade,  which,  according 
to  Kelvin,  cannot  simply  end  along  with  the  blade  at  its  tip.  At  the  blade  tip,  the 
circulation  continues  into  a  vortex,  which  is  bent  and  dragged  behind  the  moving  blade. 
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Figure  1.1:  Tip  Leakage  Flow  and  Vortex  Illustration 


1.1.2:  Detriments  to  Performance 

Tip  leakage  flows  and  their  resulting  vortices  are  detrimental  to  turbomachinery 
performance  for  multiple  reasons: 

1 )  Decay:  The  vortex  will  ultimately  decay  away,  primarily  by  turbulent  dissipation.  The 
energy  of  the  vortex  is  converted  into  internal  energy  and  cannot  be  recovered. 

2)  Underturning:  Flow  turning  is  the  primary  means  of  transferring  mechanical  energy  in 
turbomachines.  The  fluid  that  is  forced  over  the  tip  of  the  blade  is  not  turned  to  the 
extent  of  the  machine's  design.  The  result  is  less  blade  loading  near  the  tip. 

3)  Blockage:  The  tip  leakage  vortex  acts  as  an  obstruction  to  the  main  flow  through  the 
blade  passage.  Therefore,  less  mass  can  flow  through  the  machine  and  less  work  is  done. 
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For  transonic  flows  through  multiple  blade  rows,  this  blockage  can  lead  to  changes  in 
the  inlet  flow  angle  of  the  blade,  as  is  described  by  Dietrichs,  et  al.  [1]. 

4)  Entrainment:  The  slower  moving  leakage  vortex  and  its  low  pressure  core  acts  to 
entrain  the  faster  passage  flow,  thus  causing  substantial  three-dimensionality  and 
deviations  from  design  conditions. 

5)  Downstream  Unsteadiness:  In  any  turbomachinery  with  multiple  stages,  a  certain 
degree  of  unsteadiness  is  always  present  within  the  flow.  Physical  features  of  a  steady 
flow  through  a  blade  row,  such  as  wakes,  vortices,  and  boundary  layers,  are  seen  as 
oscillations  in  the  relative  flow  by  the  row  immediately  downstream.  This  is  due  to  the 
relative  motion  between  all  adjacent  blade  rows.  The  resulting  fluctuation  in  the  pressure 
and  velocity  fields  result  in  two  major  performance  detriments  that  are  of  design  concern. 
First,  the  unsteady  flow  field  means  that  blade  loading  will  not  be  steady,  but  rather 
oscillating  slightly  at  the  upstream  row’s  passing  frequency.  Although  the  magnitude  of 
this  fluctuation  is  small,  it  is  still  of  concern  from  a  fatigue  standpoint.  It  is  most  often 
desirable  to  design  turbomachines  as  lightweight  as  possible,  yet  still  able  to  endure  these 
fatiguing  stresses.  This  is  particularly  true  in  military  aircraft  applications,  which  tend  to 
call  for  operation  at  high  loading  over  long  periods  of  time. 

Second,  oscillating  pressure  due  to  upstream  relative  motion  is  the  largest  cause  of  noise 
within  turbomachinery.  This  is  particularly  true  in  compressors.  Here,  the  relative 
motion  between  the  rotor  and  the  casing  enhances  tip  vortices.  Also,  compressor 
boundary  layers  and,  thus,  wakes,  tend  to  be  larger  because  of  the  adverse  pressure 
gradients.  Both  of  these  effects  result  in  larger  pressure  oscillations  and  more  noise. 
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Noise  is  often  of  concern  to  designers  as  a  comfort  issue,  as  well  as  in  situations  where 
public  health  codes  require  limits  on  noise  emissions. 

6)  Complex  Heat  Transfer:  For  turbomachines  that  operate  at  high  temperatures,  such  as 
aircraft  engine  turbines,  the  tip  leakage  flow  increases  the  complexity  of  the  heat  transfer 
to  the  blade  and  outer  casing.  Although  heat  transfer  aspects  are  beyond  the  scope  of  this 
study,  they  are  shown  here  to  be  a  great  challenge  that  high  temperature  turbine  designers 
must  overcome.  In  such  a  turbine,  the  tip  leakage  vortex  rolls  up  along  the  suction 
surface  of  the  blade  near  to  the  tip.  This  acts  to  entrain  hot,  freestream  fluid  in  toward  the 
suction  surface  of  the  blade.  Also,  the  tip  leakage  flow  will  expose  the  top  surface  of  the 
blade  to  high  temperatures.  Finally  the  impingement  of  the  leakage  vortex  acts  to  break 
up  film  cooling  flows  and  bring  fi'eestream  gas  in  contact  with  blade  surfaces.  This  is 
true  for  both  the  blade  row  responsible  for  the  vortex  and  the  downstream  row,  where  the 
relative  xmsteadiness  of  the  situation  only  enhances  the  complexity  of  the  mixing. 

1.2:  Background 

The  term  “Tip  Desensitization”  has  been  defined  by  industry  and  academia  to 
represent  efforts  aimed  at  reducing  or  eliminating  the  leakage  flow  and  vortex  in  axial 
flow  turbines.  This  section  gives  an  overview  of  past  tip  desensitization  efforts  and 
explains  the  ciurent  understanding  of  the  physics  of  the  flow  field  inside  the  turbine  rotor. 

1.2.1:  Past  Tip  Desensitization  Efforts 

Although  there  have  been  a  niunber  of  past  studies  dealing  with  tip 
desensitization,  the  final  results  presented  tend  to  be  inconclusive.  A  few  suggested 
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concepts  in  this  area  re  the  tip  winglets,  tip  aerodynamic  modification  (blade  profile 
change)  for  reducing  leakage,  non-uniform  tip-gap  (adverse  pressure  gradient)  to  diffuse 
the  leakage  flow  structure,  normal/angled  injection  into  squealer  tip  regions,  and 
increased  turbulence  intensity  upstream  of  the  rotor  for  weakening  the  tip  vortex. 

Some  of  the  earlier  efforts  to  reduce  the  tip  clearance  effects  are  summarized 
below.  It  should  be  state  here  that  none  of  these  devices  have  been  tested  in  the  hostile 
environment  of  a  rotor  or  with  relative  motion  between  the  blade  and  the  wall. 

1.  Some  of  the  early  attempts  to  reduce  the  tip  leakage  flow  through  slots  in  the 
blade  tip  and  flow  fence  were  unsuccessful  (Lakshminarayana  and  Horlock,  [2]). 

2.  The  pressure  side  winglet  provided  improved  performance  in  the  simulated  flow 
in  a  channel  (Wadia  and  Booth,  [3]).  This  contradict  the  recent  investigations  in  a 
cascade  by  Yaras  and  Sjolander  [4],  who  concluded  that  the  pressure  side  winglet 
was  least  effective  and  the  double  sided  winglet  was  most  effective. 

3.  Radiusing  and  contouring  the  leakage  gap  geometry  prevented  the  formation  of  a 
separation  bubble  and  decreased  the  losses  (Bindon  and  Morphis,  [5]),  but  the 
resulting  thin  tip  section  may  not  be  suitable  for  a  high  temperature  application. 

4.  The  numerical  simulation  by  Chen  et  al.  [6]  indicates  that  the  two-squealer  tip 
combined  with  the  jet  injection  at  the  correct  angle  to  the  leakage  flow  reduces  the 
losses.  They  also  indicate  that  the  losses  are  very  sensitive  to  the  geometry  and 
location  of  blowing. 
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5.  Axially  non-uniform  tip  clearance  in  a  compressor  improved  the  losses  (Mohan 
and  Guruprasad,  [7]).  It  is  not  known  whether  this  concept  is  valid  for  turbine 
flow. 

6.  A  recent  investigation  by  Staubach  et  al.  [8]  indicates  that  desenstization  can  be 
achieved  through  proper  design  of  the  upstream  stator  and  regulation  of  the 
spanwise  flow  distribution  in  the  rotor  blade. 

7.  A  very  thorough  overview  of  tip  leakage  research  in  turbines,  both  experimental 
and  computational,  can  be  found  in  the  lecture  series  by  Sjolander  [9]. 

1.2.2:  Ciurent  Understanding 

A  comprehensive  numerical  study  by  Tallman  and  Lakshminarayana  [10,11] 
revealed  new  insight  into  the  physics  of  the  leakage  flow,  leakage  vortex,  and  additional 
secondary  flows  present  in  these  turbines.  Some  of  die  major  conclusions  reached  from 
this  study  are  outlined  below,  and  labeled  in  Fig.  1.2. 

1 .  For  tip  clearance  heights  representative  of  a  real  turbine,  the  leakage  vortex  (Fig. 
1.2  Lab.  i)  is  not  caused  by  a  leakage  jet  shearing  across  the  passage  flow,  as  was 
previously  believed. 

2.  Most  of  the  leakage  flow  through  the  gap  region  does  not  roll  into  the  leakage 
vortex  (Fig.  1.2  Lab.  ii).  Only  the  fluid  that  passes  through  the  gap  very  near  to 
the  blade  tip  ends  up  with  the  leakage  vortex  (Fig.  1.2  Lab.  iii).  The  rest  of  the 
leakage  flow  joins  with  the  passage  flow,  but  still  contributes  to  additional 
secondary  flow  effects  inside  the  passage  (Fig.  1.2  Lab.  iv). 
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3.  Vorticity  convection  by  the  leakage  flow  from  the  gap  and  into  the  passage  plays 
an  important  roll  in  the  physics  of  the  secondary  flow  field.  Fluid  passing  through 
the  gap  near  the  blade  tip  exits  the  gap  with  rotation  in  a  direction  closely 
corresponding  to  the  rotation  in  the  leakage  vortex. 

z 


Figure  1.2:  Secondary  Flows  in  an  Axial  Flow  Turbine 

4.  A  strong  secondary  flow  toward  the  suction  side  of  the  passage  is  present  near  the 
outer  casing  in  the  turbine  (Fig.  1.2  Lab.  v).  This  secondary  flow  resembles  a 
passage  vortex  near  the  casing,  but  is  primarily  caused  by  the  relative  motion  of 
the  outer  casing  with  respect  to  the  blade.  This  secondary  flow  acts  to  obstruct 
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the  leakage  flow  throu^  the  gap,  particularly  near  the  leading  edge  of  the  gap, 
and  acts  to  reduce  the  leakage  flow  and  vortex.  This  secondary  flow  was  given 
the  name  “Near-Casing  Secondary  Flow”  by  Tallman  and  Lakshminarayana 
[10,1 1],  and  is  also  adapted  here. 

5.  Two  additional  secondary  flows  were  also  defined  by  Tallman  and 
Lakshminarayana  [10,11].  The  “Blockage  Secondary  Flow”  results  as  the 
passage  flow  is  obstracted  around  the  leakage  vortex  and  the  near-casing 
secondary  flow  (Fig.  1.2  Lab.  vi).  The  “Gap  Entrainment  Secondary  Flow”  is 
caused  by  the  entrainment  of  the  leakage  fluid  into  the  gap  along  the  pressure 
surface  side  (Fig.  1.2  Lab.  vii). 

1.3:  Objectives  and  Direction  of  the  Research 

The  current  research  being  presented  here  is  a  continuation  of  the  original 
Tallman  and  Lakshminarayana  study  [10,1 1].  The  objective  of  the  current  research  is  to 
study  the  effect  of  modified  blade  tip  designs  on  the  leakage  vortex  for  axial  flow 
turbines.  Computational  Fluid  Dynamics,  or  CFD,  is  used  as  the  tool  for  study  and 
analysis  of  the  flow  field.  Ultimately,  these  modified  blade  tip  designs  are  aimed  at 
reducing  the  leakage  vortex  and  its  negative  effects. 

1.3.1:  Blade  Tip  Designs  Simulated 

Through  some  creative  grid  stretching  in  the  spanwise  direction,  various  degrees 
of  tip  chamfering  were  added  to  the  turbine  blade  geometry  simulated  in  the  original 
study.  For  clarity,  a  chamfered  blade  tip  is  illustrated  in  Fig.  1.3.  These  chamfered  blade 
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tip  designs  take  advantage  of  the  flow  physics  insight  gained  in  the  original  Tallman  and 
Lakshminarayana  study  [10,11],  and  are  aimed  at  reducing  the  magnitude  of  the  leakage 
vortex.  Three  cases  are  simulated  for  comparison  with  a  regular,  horizontal  tip  baseline 
case.  The  baseline  case  is  taken  from  Tallman  and  Lakshminarayana  [11].  The  three 
modified  cases  and  their  design  reasoning  are  outlined  below. 

1.  T.Radinp  KHpe  Chamfer:  The  near-casing  secondary  flow  is  seen  in  Fig.  1.2  to  act 
as  an  obstruction  to  the  leakage  flow.  By  increasing  the  gap  height  on  the 
suction  side  near  the  leading  edge,  the  near-casing  secondary  flow  may  enter  the 
gap  region  and  fiirther  inhibit  the  leakage  flow. 

2.  Trailing  Fdpe  Chamfer:  By  increasing  the  gap  height  on  the  suction  side  from  the 
mid-chord  to  trailing  edge  of  the  blade,  the  increasing  flow  area  in  the  gap  may 
act  as  a  diffuser.  This  could  reduce  the  leakage  flow’s  velocity  and,  thus,  it’s 
influence  on  the  leakage  vortex  and  secondary  flow  field  outside  the  gap. 

3.  Full  Chamfer:  Both  of  the  first  two  cases  applied  together. 

In  each  of  the  three  cases,  the  tip  clearance  height  along  the  pressure  side  of  the 
blade  is  held  at  a  uniform  1 .0%  of  the  blade  span  (1.86  mm).  Therefore,  the  minimum  tip 
clearance  is  the  same  for  all  of  the  cases.  Also  for  each  of  the  three  cases,  the  entire 
blade  tip  thickness  is  not  chamfered.  On  the  pressure  side  of  the  gap  region,  the  tip 
remains  horizontal  for  approximately  10%  of  the  distance  across  the  gap.  This  short 


horizontal  inlet  region  is  illustrated  in  Fig.  1.3.  It  is  included  to  avoid  sharp  comers  on 
the  blade  tip,  which  could  cause  flow  separation  and  would  be  difficult  to  cool. 
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Figure  1.3:  Chamfered  Blade  Tip  Illustration 


For  all  three  cases,  the  angle  of  the  chamfer  angle  is  7.7  degrees  from  horizontal. 
This  angle  corresponds  to  a  suction  side  clearance  height  of  4%  blade  span  (7 .44  mm)  at 
the  blade’s  thickest  location.  Note  that  the  clearance  height  will  vary  along  the  suction 
surface  due  to  the  blade’s  varying  thickness,  but  remains  constant  (1.86  mm)  along  the 
blade’s  pressure  surface.  In  each  of  the  three  cases,  the  chamfering  angle  is  zero  degrees 
at  the  inlet  and  exit  of  the  computational  domain.  Near  the  leading  edge  of  the  blade,  the 
chamfer  angle  is  linearly  increased  from  zero  to  7.7  degrees.  Near  the  trailing  edge  of  the 
blade,  the  chamfer  angle  is  linearly  decreased  back  from  7.7  to  zero  degrees.  For  the  full 
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chamfer  case,  the  chamfer  increase  takes  place  upstream  of  the  blade  leading  edge  and 
the  chamfer  decrease  takes  place  downstream  of  the  blade  trailing  edge.  For  the  leading 
edge  chamfer  case,  the  chamfer  increase  takes  place  upstream  of  the  blade  leading  edge 
and  the  chamfer  decrease  takes  place  on  the  blade  between  47.4%  and  63.4%  axial  chord. 
For  the  trailing  edge  chamfer  case,  the  chamfer  increase  takes  place  between  38.1%  and 
70.3%  axial  chord  and  the  chamfer  decrease  takes  place  downstream  of  the  blade  trailing 
edge.  The  three  chamfered  tips  are  illustrated  in  chapter  3  (Fig.  3.1). 

1.3.2:  Reasons  for  Nmnerical  Simulation 

Computational  Fluid  Dynamics  (CFD)  was  chosen  as  the  means  of  investigation 
for  this  study.  CFD  is  logistically  better  suited  to  the  study  of  the  fine  details  of  the  tip 
leakage  flow  and  vortex  than  experimental  measurement.  Even  with  present  advanced 
measurement  techniques,  such  as  laser  Doppler  velicometry,  measuring  the  flow  inside 
the  tip  clearance  gap  remains  an  impossibility.  CFD  offers  a  tool  for  looking  at  the  flow 
inside  areas  too  small  for  measurement.  Also,  CFD  is  considerably  less  expensive  for 
testing  multiple  prototype  blades  than  manufacturing  test  blades  for  experiment.  The 
numerical  procedure  that  was  used  in  this  study  is  outlined  and  validated  in  chapter  2,  and 
has  been  previously  validated  by  Basson  and  Lakshminarayana  [13]. 

1.3.3:  Choice  of  Blade  Geometry 

A  linear  turbine  cascade  blade  originally  studied  by  Bindon  [5,12]  was  chosen  as 
the  geometry  to  be  used  in  this  study.  This  blade  is  shown  in  Fig.  1.4,  and  an  overview  of 
its  geometry  is  given  in  Table  1.1.  This  is  the  same  turbine  cascade  that  was  used  in  the 
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original  Tallman  and  Lakshminarayana  study  [10,11].  The  cascade’s  profile  does  not 
change  in  the  spanwise  direction,  isolating  the  leakage  flow  and  vortex  from  the  effects 
of  sweep,  lean,  or  any  other  complicated  geometrical  parameters. 

Two  modifications  have  been  made  to  the  linear  turbine  cascade  of  Bindon  [5,12]. 
First,  the  tip  clearance  height  is  reduced  from  2.5%  to  1.0%  of  the  blade  span.  Second, 
unlike  a  cascade,  the  effect  of  casing  relative  motion  is  included  in  the  current 
simulations,  using  a  wall  velocity  of  31.0  m/s  toward  the  suction  side  of  the  passage. 
This  wall  velocity  corresponds  to  a  flow  coefficient  value  of  0.45,  where  the  flow 
coefficient  is  defined  as  the  inlet  axial  velocity  divided  by  the  moving  blade’s  velocity. 
These  values  of  tip  clearance  height  and  flow  coefficient  are  more  typical  of  a  low- 
pressure  turbine. 


Table  1.1:  Turbine  Case  Information 


Chord  Length  =  Span  Length 

0.186  m 

Pitch  Length 

0.1323  m 

Inlet  Camber  Angle 

0.0° 

Exit  Camber  Angle 

68.0° 

Tip  Gap  Height 

1.0%  Chord 

Inlet  Flow  Velocity 

13.82  m/s 

Inlet  Boimdary  Layer  Thickness 

10%  Span 

Inlet  Turbulence  Intensity 

Approx.  2%  Inlet  Velocity 

Outer  Casing  Relative  Velocity 

-31.0  m/s  tangential 

Working  Fluid 

Air 
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Axial  Direction  /  Chord  Length 


Figure  1.4:  Blade  Profile  and  Velocities  Used 


1.4:  Report  Overview 

Chapter  1  of  this  report  has  introduced  the  leakage  vortex  concept,  its  ill  effects, 
our  current  imderstanding  of  the  leakage  vortex,  and  the  current  research  program. 
Chapter  2  presents  the  numerical  method  that  is  used  for  the  four  different  simulations  in 
this  report.  In  addition,  some  code  validation  results  from  a  separate  turbine  cascade 
simulation  are  presented  in  chapter  2,  which  show  that  the  numerical  method  is  capable 
of  predicting  complicated  turbine  secondary  flows.  Chapter  3  presents  the  simulation 
results  for  the  three  modified  blade  tips  discussed  in  section  1.3.  For  each  of  the  three 
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cases,  results  are  compared  directly  with  those  from  the  baseline,  horizontal  blade  tip 
case.  Chapter  4  presents  the  conclusions  found  in  this  research,  as  well  as  ideas  for  future 
studies  that  build  on  the  current  research. 

Some  common  terms  that  are  used  throughout  the  report  are  defined  here  for 
clarity.  The  axial  (x),  tangential  (y),  and  spanwise  (z)  directions  are  referred  to 
throughout  this  report  according  to  their  common  definitions.  These  directions  are 
illustrated  in  Fig.  1.3.  The  term  “2D  Streamwise”  is  used  to  refer  to  the  direction  of  the 
2D,  inviscid,  potential  flow  through  the  passage.  The  term  “Camber  Direction”  is  used  to 
refer  to  the  local  direction  of  the  blade’s  camber  line.  The  term  “Secondary  Flow”  is 
used  generically  to  represent  any  flow  deviation  from  the  2D  streamwise  direction,  and 
should  not  be  confused  with  the  term  “Passage  Vortex”. 


CHAPTER  2:  NUMERICAL  METHOD  REVIEW  AND  VALIDATION 


The  full  Reynolds  averaged  Navier-Stokes  equations  are  employed  in  this 
research  to  model  the  flow  through  a  turbine  cascade.  The  computer  code  used  in  this 
study  uses  a  pressure-correction  method  for  solving  conservation  equations  for  mass, 
momentum,  and  turbulence  quantities.  A  pressure  correction  equation  is  used  to  correct 
the  pressure  field  such  that  the  continuity  equation  is  also  satisfied.  Because  the  method 
uses  a  non-staggered  grid,  additional  artificial  dissipation  terms  are  added  to  these 
equations  to  aid  in  smoothing  out  errors  and  to  prevent  numerical  oscillations  of  the 
centrally  differenced  pressure  equation. 

Chapter  2  is  organized  as  follows:  Section  2.1  presents  the  governing  equations  in 
both  their  standard  and  discretized  forms.  The  complete  derivation  of  these  equations  can 
be  found  in  the  thesis  written  by  Basson  [14],  and  a  more  general  text  on  the  method  is 
that  of  Patankar  [15].  Section  2.2  presents  the  grid  that  was  used  and  compares  the  grid 
to  its  predecessors.  Section  2.3  compares  the  results  of  an  initial  computation  with 
experimental  data  for  validation  of  the  code. 

2.1:  Governing  Equations 

The  Reynolds  averaged  equations  of  mass  and  momentiun  conservation  are  given 
in  equations  2.1  and  2.2,  respectively.  Standard  tensor  notation  is  used  for  their 
presentation,  where  i  represents  1, 2,  or  3;  j  represents  smnmation  over  1,  2,  and  3;  and  1, 
2,  and  3  correspond  to  the  x,  y,  and  z  directions,  respectively.  Steady  (in  the  mean), 
constant  density  flow  conditions  are  assumed  in  their  derivations. 
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d  Ui  d  UfUj 

dt  dxi 


I  dp  d  [  5  m,.  I  _ — tt 

— - u  — -  +  — -  +pu]Ui 

p  dXj  dxj  dxj  J 


Tiirbulence  closure  is  obtained  via  a  standard  k-e  model,  where  the  Reynolds 
stress  term  at  the  end  of  equation  (2.2)  is  modeled  as: 


Conservation  equations  for  k  and  s  therefore  must  also  be  solved.  The  k  and  s 

equations  are  given  in  equations  (2.4)  and  (2.5),  respectively. 

dpk  dupk_  (^•^) 

dt  dxi 

d  {(  ^dk^  duifdui  5m.")  2pk 

dx,)  dx\dxj  dxj  (yj 


ds  5m,£- 

—  +  — t— 
dt  dXj 


//,  ^  sduJdUi  duj 

•—  -  +Qp, - ^  — -  +  ^r^ 

kdxj^dxj  dxi 


k 


The  near  wall  damping  functions,  f2,  f^,  and  fe  are  from  the  k-s  model  by  Chien 
(1982).  The  functions  and  empirical  constants  used  in  the  model  are  given  below. 

4-1  ^  4  1  0.8"  {yJ 


C^=0.09  cr,  =1.3  C,  =1.35  ^=1.8  Cj  =0.0115  C4=0.5 
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The  continuity,  momentum,  and  turbulence  equations  can  all  be  represented  in 
the  following  general  convection-diffusion  form: 

(2.6) 


d  p(t>  g  puj<i>  _  d 
dt  ^ 


dxj  dxj 


- 


dxjj 


Here,  ^  represents  1,  Ui,  k,  or  s  for  equations  (2.1),  (2.2),  (2.4),  and  (2.5),  respectively. 
represents  the  corresponding  diffusion  coefficient,  which  can  be  laminar  or  turbulent.  S,j, 
is  the  source  term. 


2.1.1 :  Computational  Domain  Transform 

For  its  solution,  equation  (2.6)  is  first  transformed  fi-om  the  real  domain  (x,y,z)  to 
the  computational  domain  (^,ti,Q  and  then  discretized  using  finite  difference 
approximations.  The  computational  domain  is  such  that  all  grid  spacings  are  unity  and 
all  grid  lines  are  orthogonal.  The  transform  is  accomplished  using  the  chain  rule  of 
calculus  and  Kramer’s  rule  of  matrix  algebra.  The  general  convection-diffusion  equation 
in  the  computational  domain  is  given  in  equation  (2.7). 


dt  07]  0(, 


d 

rO 

(  d<i>  ^  .  dA 

d 

+  - 

rO 

[  d4  'dfj 

df] 

^  dt]  d^)_ 

+ 


_a_ 


dt] 


+  JS^ 


(2.7) 
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The  new  geometrical  terms  in  equation  (2.7)  are  defined  as  follows: 


Gi  =Uj- 


Gj  =  Uj 


'  "  dx, 


KK.  0=^^  y=KK.  0  0  0 

dxj  dxj  ^  dxj  dxj  ^  dxj  dxj  ‘  dxj  dxj  ‘  dxj  dxj  *  dxj  dxj 


dxJdx^  dx,  dxi  f  dx^  dx^  8x2  ^^^3^  ^  dx^  f  dx^  dx^  dx^  dx^'^ 

'^ydT]  dC  dC  di]J  d7i[d^  dC  j  dc[d^  drj  dtj  j 


The  terms  Gi,  G2,  and  G3  are  the  contravariant  velocities  in  the  r|,  and  C,  directions, 
respectively.  J  is  the  Jacobian  of  the  coordinate  transform,  and  a,  P,  y,  0i,  82,  and  83  are 
geometric  parameters.  The  new  geometrical  terms  are  calculated  using  finite  difference 
approximations  from  the  grid.  All  differences  in  y\,  and  ^  are  unity. 

In  addition,  two  artificial  terms  are  added  directly  to  equation  (2.7)  for  the 
momentum  and  turbulence  equations.  In  the  transformed  coordinates,  these  terms  are: 


These  terms  act  to  artificially  diffuse  errors  out  of  the  solution.  The  coefficients  82  and  84 
are  referred  to  as  the  2"^  and  4*  order  artificial  viscosity,  respectively,  and  are  specified 
as  input  to  the  code.  In  any  final  solution,  82  and  84  should  be  set  as  low  as  possible  to 
avoid  corrupting  the  solution  with  these  nonphysical  diffusion  effects. 
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2.1.2:  Finite  Difference  Form 

A  control  volume  around  a  grid  point  P  is  used  for  the  discretization  of  equation 
(2.7).  Such  a  control  volume  is  shown  in  figure  2. 1 . 


T 


Upper  case  letters  represent  grid  points  in  the  computational  domain  (N  =  north),  while 
lower  case  letters  represent  the  faces  of  the  control  volume.  The  final  form  of  the 
discretized  general  convection-diffusion  equation  will  be  as  follows: 


The  subscripts  indicate  the  grid  point  location. 
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For  stability  reasons,  it  is  beneficial  if  the  above  equation  can  be  solved  in  an 
implicit  manner.  The  SOURCE  term  and  “a”  coefficients  are  always  found  explicitly 
from  grid  geometry  information  and  current  values  of  (j).  New  values  of  (|)  are  found 
using  a  tridiagonal  matrix  solver  in  a  semi-implicit  manner.  For  example,  the  terms  aN(|)Nj 
as(|)s,  aT<l)T,  and  aB<|)B  will  be  foimd  explicitly  and  subtracted  from  the  SOURCE  term  for 
all  grid  points  along  a  line  of  constant  4-  Since  <|)p  is  equivalent  to  (|)e  for  the  grid  point  to 
the  left  and  (|)w  for  the  grid  point  on  the  right,  including  boundary  conditions  will  result  in 
a  set  of  n  equations  and  n  unknowns.  These  equations  make  up  an  easily  solvable 
tridiagonal  matrix,  with  n  values  of  (|)p  as  the  unknowns.  After  solving  for  new  values  of 
(t»p,  the  process  is  repeated  along  lines  of  constant  11  (<t)p,  <|)n,  and  (|)s  treaded  implicitly), 
and  constant  ^  (<|)p,  (j)T,  and  (|)b  treaded  implicitly), 

2"**  order  accurate  central  differencing  is  used  for  both  convection  and  real 
diffusion  terms.  For  example,  the  convection  derivative  in  the  ^  direction  is  discretized 
as  follows: 

-^\pJG,(l>]=  A^  =  A7  =  A^  =  1 

The  values  at  the  control  volume  faces  come  from  linear  interpolation.  For  example: 
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The  diffusion  terms  are  modeled  similarly,  as  follows: 


In  equation  (2.10),  derivatives  at  the  control  volume  faces  are  central  difference 
approximated.  Cross  derivative  terms  are  calculated  on  either  side  of  the  face  and  then 
averaged.  Examples  of  these  approximations  are  given  below. 

5^''  _  i^E  ~^p)  ~^>r) 

_  1  {^N  i^NE  ~  ^SE  ) 

2(Av)t  2(A^)  2(A7)  J 

_  1  ~^s)  ^  ^NW  ~  ^SW  y\ 

■  2(A7)t  2(A7)  2(A7)  J 

It  can  be  seen  here  that  the  modeling  of  cross  derivative  diffusion  terms  will 
require  additional  grid  points  other  than  the  ones  shown  in  figure  2.1.  In  order  to  solve 
the  final  equations  with  a  tridiagonal  matrix  solver,  these  points  will  need  to  be  treated 
explicitly. 

The  2"^  order  artificial  diffusion  term  is  modeled  in  the  same  manner  as  the  real 
diffusion  terms.  The  4^'’  order  artificial  term  is  calculated  by  the  central  difference 
approximation  of  two  3*^  order  artificial  derivatives,  as  is  shown  in  equation  2.1 1. 
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—  ^WW  ^  ^W  '^^P  ■*■  ^£ 


|/7JG,| 


^ay^ 


/?/ Gi 


yy^ 


^ay^ 


=  |p/Gi|^(-^,y  +3^^;,  -34  +^££)-|/?/ -3^i>  +<I>e) 


No  cross  derivatives  are  present  in  the  artificial  terms.  The  third  order  derivative 
approximations  in  equation  (2.1 1)  can  be  written  as: 


(2.12) 


ay 

a#^ 


=  -(!>„  +'i<j)p  -3^£  =  ir^W  ~^P  +^£  +^^£e)“4(“^P  '^^e) 


Again,  the  use  of  a  tridiagonal  solution  technique  requires  that  the  terms  in  the  first 
brackets  on  the  right  side  of  equation  (2.12)  must  be  treated  explicitly. 

The  final  discretized  form  of  the  general  convection-diffusion  equation  is 
presented  in  equation  (2.13).  S,|,  is  the  source  term  for  the  particular  equation  being 
solved:  either  the  pressure  term  in  the  momentum  equations  or  the  production  and 
dissipation  terms  in  the  turbulence  equations.  These  terms  are  first  transformed  to  the 
computational  domain  by  the  method  described  in  section  2.1.1,  and  then  discretized 
using  central  differencing.  Their  treatment  is  always  explicit. 
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♦  J 

Up  (j)^  +  Cljy  +  Clff  +  fl J  <j>j  +  Ug  +  S^ 

where 

a,  =  -iO’/G,),  +(rVa),  +ife  +  «.)|(p/G,),| 
a,  =  +(rVa).  +ife  +^,)|(p/G,),| 

a„  =  -i(p/G,),  +(rVy9),  +ife  +i.)|(p/G,).| 

a,  =  i(p/G,),  +(rv4  +|fe +«4)|(/’-^GJ,| 
ar  =-|(p/G,),  +{r'>/r),  +ife  +^4)|(/’-/G>),| 
a,  =  -jC^/G,),  +(T*Jr\  +ife  +«4)|(/’-/G,).| 

flp  =  +  At-  + 

+  [(/?/  (j,)g  +(/?  JG2)„  -(/?/(j2)j  “(P'^^s)*] 

5*  =  JS^  +  [Cro5J  Derivative  Terms\+  [Explicit  Artificial  Viscosity  Terms] 


2.1.3:  Pressure  Correction  Equation 

The  continuity  equation  does  not  lend  itself  well  to  solution  techniques  in 
conjunction  with  the  momentum.  The  reason  for  this  is  the  lack  of  a  pressure  term  in  the 
continuity  equation.  An  alternative  approach  is  to  solve  a  pressure  or  pressure  correction 
equation  derived  directly  from  the  discretized  continuity  and  momentum  equations.  This 
is  the  essence  of  the  Semi-Implicit  Method  for  Pressure  Linked  Equations  (SIMPLE) 
scheme  developed  by  Patankar  [15]. 
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The  pressure  correction  equation  can  be  derived  as  follows:  First,  the 
discretized  momentum  equations  are  substituted  into  the  discretized  continuity  equation, 
resulting  in  an  expression  for  pressure  in  terms  of  velocity  at  any  of  the  grid  points. 
Then,  the  exact  pressure  in  this  equation  is  replaced  by  the  sum  of  the  current  guessed 
pressure  value,  p*,  and  a  pressme  correction  value  p^.  If  the  current  velocity  and  pressure 
values  satisfy  the  momentum  equations,  then  these  equations,  using  p*  and  not  p,  can  be 
subtracted  from  the  pressure  correction  equation.  The  final  pressure  correction  equation 
is  given  in  equation  (2.14). 


^pP'p  =  hew  1  P'e  +  hew  L  P'w  +  hns  I  P'n  +  hns  1  P's  +  ho,  1  P't  +  htb  L  P'b 
-[pJG\+[pJG’X-[pJG’X+[pJG'X-[pJG[\  +[pJG'X 


EAe,—+e^— 


E^J\ 


E,J\ 


0. 


XeX 


E^J\ 


94  dC 


Jt 


E^J\ 


dtt 


where 


E  = 


^P  ^nb^nb 


,  nb  =  E,W,N,S,T,B  and  a  =  coef.jrom  (2.13) 


a^  =  pEJ^a  a„^=pEJ^p  a^=pEJ^Y  E^=pEJ 

V/>=^|2««6(v„6-Vp)+5, 


(2.14) 
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Uf  =  correct  =  u*  +  u^  =  current  +  correction  (i  =  1,2,3) 
M;  =  correct  =  u]  +  =  current  +  correction  (/  =  1,2,3) 


=  u  .  — ^ 

*  c 

*  drj 

-Ui — - 

Gl- 

=  u. 

‘  dXj 

^  dxj 

li 

G'.  -- 

-u] 

A 

G's  = 

^  dxj 

2 

^  dxj 

In  addition,  a  4*  order  artificial  dissipation  similar  to  that  in  equation  (2.8)  is 
included  in  the  pressure  correction  equation.  This  is  necessary  because  the  central 
differencing  in  the  pressure  correction  equation  leads  to  an  odd-even  decoupling  effect, 
which  is  described  in  detail  in  the  text  by  Patankar  [15].  The  coefficient  8pw  is  specified 


directly  and  should  be  chosen  as  small  as  possible.  The  additional  term  is  shown  in 
equation  (2.15),  and  its  discretization  is  similar  to  that  of  the  4*  order  artificial  term  in 
the  general  convection-diffusion  equation. 


1  \  ^  (  Af  ^ 

drf  j 


(2.15) 


2.1.4:  Boxmdarv  Conditions 


There  are  four  types  of  boimdaries  to  the  computation  region.  The  boundary 
conditions  that  were  implemented  on  these  four  types  of  boundaries  are  explained  below. 
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Inlet: 

Inlet  velocity  and  turbulence  profiles  are  read  directly  into  the  code.  For  this 
thesis,  the  origin  of  these  profiles  is  dependent  upon  the  case  being  run.  For  the  code 
validation  case  in  chapter  two,  the  boundary  layer  thickness  and  fireestream  turbulence 
intensity  were  chosen  to  agree  with  the  experimental  results  being  compared  with.  These 
values  were  approximately  10%  of  the  chord  length  for  the  boundary  layer  tihickness  and 
2%  of  the  inlet  velocity  for  turbulence  intensity.  The  inlet  turbulence  length  scale  was 
unknown  from  the  experiment  and  chosen  as  approximately  10%  of  the  chord  length. 
For  the  boundary  layer  profile,  the  commercial  CFD  code  Fluent  was  used  to  compute  a 
2D  flat  plate  boimdary  layer  with  the  same  ffeestream  conditions  as  are  given  above.  The 
profile  at  the  proper  thickness  was  then  used  for  the  inlet  conditions.  For  the  boimdary 
layer  computation,  a  full  Reynolds  stress  turbulence  model  was  used,  which  was  most 
likely  calibrated  using  experimental  flat  plate  boundary  layer  measurement  data. 

The  Fluent-calculated  inlet  profile  was  also  used  for  the  reduced  clearance  case 
presented  in  chapter  3.  In  the  cases  in  chapter  4,  a  more  realistic  inlet  profile  was  used. 
This  second  profile  was  computed  in  the  manner  explained  in  section  1.3. 

Pressure  is  assumed  to  be  unchanging  in  the  flow  direction  at  the  inlet.  Therefore, 
the  inlet  should  be  located  far  enough  upstream  of  the  blades  to  ensure  that  this 
assumption  is  valid.  For  this  study,  the  inlet  was  placed  0.6  chord  lengths  upstream  of 
the  blade  leading  edge. 
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Exit: 

All  convection  and  diffusion  terms  for  velocities  and  turbulence  quantities  are 
computed  directly  at  the  exit  plane  with  no  direct  enforcement  of  a  boundary  condition. 
Static  pressure  at  the  exit  is  directly  specified  as  zero.  To  ensure  that  this  assumption  is 
valid,  the  exit  plane  must  be  chosen  far  downstream  of  any  blading.  For  this  study,  the 
exit  plane  is  located  approximately  two  chord  lengths  downstream  of  the  blade  trailing 
edge. 

Solid  Surfaces: 

Setting  velocities  and  the  turbulent  kinetic  energy  to  zero  on  all  solid  boundaries 
enforces  the  no-slip  condition.  The  exception  to  this  is  at  the  outer  casing,  where  the 
velocity  in  the  tangential  direction  is  set  to  a  constant,  in  order  to  simulate  endwall 
relative  motion  in  the  planar  geometry.  For  relative  endwall  motion  in  an  annular 
geometry,  two  Cartesian  velocity  components  must  be  specified,  with  their  values 
depending  on  the  angular  location  on  the  wall.  The  near-wall  damping  terms  of  Chien 
also  require  that  the  dissipation  rate,  s,  be  set  to  zero  on  all  solid  boimdaries. 

Wall  pressure  is  determined  by  setting  the  pressure  gradient  normal  to  the  surface 
equal  to  zero.  In  order  to  ensure  the  validity  of  this  assumption,  as  well  as  that  of  the 
turbulence  model,  the  first  grid  point  out  firom  each  solid  boundary  is  placed  at  a  y^ 


distance  of  order  one. 
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Periodic  Boundaries: 

Outside  of  extreme  operating  conditions,  the  flow  through  a  turbomachinery  blade 
row  is  periodic,  that  is,  the  flow  at  any  location  in  the  blade  row  is  identical  to  the  flow  at 
a  point  N  blade  spacings  away  in  the  tangential  direction,  where  N  is  some  integer.  This 
simplifies  turbomachinery  computations  considerably,  since  the  flow  only  needs  to  be 
solved  for  across  one  blade  spacing.  The  periodic  definition  stated  above  is  enforced 
numerically  for  all  boundary  points  that  are  not  solid,  inlet,  or  exit.  Annular  geometry 
requires  that  cylindrical  velocities  be  made  equivalent  and  not  Cartesian  velocities. 

2.2:  Grid 

The  CFD  code  used  in  this  thesis  was  specifically  written  for  H  type  grids  with 
additional  grid  points  embedded  into  the  tip  gap.  This  type  of  grid  is  common  in 
turbomachinery  computation  where  the  use  of  a  multiblock  grid  is  not  possible.  One 
problem  that  is  common  with  embedded  H  type  grids  is  the  excessive  number  of  grid 
points  present  within  a  small  location  upstream  and  downstream  of  the  leading  and 
trailing  edges,  respectively.  Additional  problems  are  grid  skewness  and  large  aspect 
ratios  in  the  gap  near  the  leading  and  trailing  edges. 

The  embedded  H  grid  used  in  this  study  is  shown  in  figures  2.1(a)  and  2.1(b). 
The  embedded  region  has  been  widened  considerably,  alleviating  the  over-clustering 
around  the  periodic  line.  The  new  grid  also  allows  for  more  efficient  clustering  of  grid 
points  along  the  blade  surfaces,  for  better  resolution  of  the  boundary  layer.  Finally,  grid 
clustering  near  the  leading  and  trailing  edges  of  the  blade  was  incorporated  in  order  to 
improve  the  average  grid  cell  aspect  ratio  in  those  regions. 
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The  biggest  drawback  to  the  modified  grid  came  in  its  generation.  An  H  grid 
developing  code  written  by  Basson  [14]  was  modified  substantially  and  then  used  to 
generate  the  region  labeled  “inner”  in  figure  2.1(a).  This  code  uses  a  combination  of 
algebraic  and  elliptic  methods,  in  conjunction  with  the  Minimal  Residual  Method,  for 
generating  smooth  grids.  The  Embedded  region,  labeled  as  such,  was  developed  using 
mostly  algebraic  techniques.  Grid  stretching  in  the  ^  direction  for  this  region  was 
accomplished  using  a  ratio  method,  which  matches  the  ratios  of  distances  between  grid 
points  along  lines  of  constant  ^  to  a  corresponding  ratio  at  the  inner  region’s  north  or 
south  boundary.  Polynomial  splines  were  used  in  this  region  in  order  to  match  the  slopes 
between  the  inner  and  embedded  regions,  ensuring  a  smooth  transition  between  grid  cells 
in  the  r\  direction.  Certain  grid  points  near  the  intersection  of  the  blade  leading  or  trailing 
edge  and  the  two  regions  were  smoothed  on  an  individual  basis. 

Because  there  is  no  variation  of  the  blade  profile  in  the  spanwise  direction,  the 
baseline  case  3D  grid  is  created  by  simply  stacking  up  the  2D  grid  in  figure  2.1.  The 
distribution  in  the  spanwise  direction  was  created  by  primarily  algebraic  methods. 
Spacing  near  the  hub,  casing,  and  tip  were  chosen  small  enough  to  give  an  approximate 
y+  value  in  single  digits.  Grid  points  were  heavily  clustered  near  the  tip,  with  28 
spanwise  grid  points  in  the  gap  and  60  spanwise  grid  points  within  six  gap  heights  of  the 
tip.  A  total  of  103  and  161  grid  points  were  used  in  the  spanwise  and  axial  directions, 
respectively.  In  the  tangential  direction,  50  and  26  grid  points  were  used  across  the  inner 
and  embedded  regions,  respectively.  Additional,  non-uniform  grid  stretching  in  the 
spanwise  direction  was  used  to  create  the  three  chamfered  grids  used  in  this  study. 
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2.2.1:  Coordinate  System 

The  Cartesian  coordinate  system  used  throughout  the  thesis  is  shown  in  figure  2.1. 
Coordinates  X,  Y,  and  Z  correspond  to  the  axial,  tangential,  and  spanwise  directions, 
respectively.  In  addition,  two  more  directions  are  referred  to  throughout  the  thesis.  The 
“2D  Streamwise”  direction  refers  to  the  direction  of  the  2D,  inviscid,  potential  flow 
through  the  passage.  The  “Principle  Normal”  direction  refers  to  the  direction  normal  to 
the  2D  streamwise  direction  in  the  blade-to-blade  plane  and  towards  the  pressure  surface. 
The  generic  term  “Secondary  Flow”  refers  to  total  flow  in  the  principle  normal  and 
spanwise  direction,  and  is  thus  the  deviation  of  the  flow  from  its  2D  designed  course. 


Figure  2.2(a):  Grid  (Every  Other  Gridline  Shown) 
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2.3:  Code  Validation 

In  order  to  validate  the  CFD  code,  and  as  a  baseline  reference,  the  flow  through 
the  turbine  cascade  described  in  table  1.1  was  computed  and  compared  with  experimental 
data  due  to  Bindon  [5,12].  Basson  [14]  and  Engrebretsen  [16]  have  carried  out  similar 
computations  with  the  same  code.  Both  studies  yielded  results  similar  to  those  found  in 
this  study.  Overall,  the  solution  method  using  k-e  turbulence  closure  is  able  to  resolve 
the  size,  magnitude,  and  location  of  the  tip  leakage  flow. 

The  pitchwise  flow  angle  is  defined  as  the  angle  between  the  total  flow  in  the 
blade-to-blade  plane  and  the  axial  direction,  and  is  positive  if  the  flow  has  a  velocity 
component  in  the  positive  tangential  direction.  The  spanwise  flow  angle  is  defined  as  the 
angle  between  the  total  velocity  in  three  dimensions  and  the  blade-to-blade  plane,  and  is 
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positive  if  the  flow  has  a  velocity  component  in  the  positive  spanwise  direction. 
Figure  2.3  compares  the  computed  and  measured  pitchwise  flow  angles.  Figure  2.4 
shows  the  same  comparison  for  the  spanwise  flow  angle.  Both  figures  are  presented  in 
three  parts:  parts  a,  b,  and  c  correspond  to  axial  locations  of  60,  80,  and  90  percent  of  the 
axial  distance  downstream  of  the  blade  leading  edge,  respectively.  For  each  part,  the 
variation  of  the  angle  in  the  spanwise  direction  is  shown  at  eight  different  pitch  locations 
ranging  from  pressure  to  suction  surfaces.  The  term  “Percent  Pitch”  is  used  to  refer  to 
the  percentage  of  the  distance  in  the  tangential  direction  from  the  pressure  surface  to  the 
suction  surface. 

At  60%  chord,  the  magnitude  of  the  leakage  vortex  location  is  found  to  be  slightly 
closer  to  the  blade  than  in  the  measured  case.  Figure  2.3(a)  shows  the  pitchwise  flow 
angles  to  be  in  good  agreement  except  at  the  80  and  90  percent  pitch  locations.  These 
locations  exhibit  the  same  trends,  but  to  a  lesser  extent  for  the  numerical  solution.  In 
figure  2.4(a),  the  same  trends  are  seen  for  the  spanwise  flow  angle.  Furthermore,  at  60, 
70,  and  80  percent  pitch  locations,  the  computed  spanwise  flow  angle  at  one  pitch 
location  appears  to  be  predicted  in  better  agreement  with  the  experimental  data  at  a  pitch 
location  fixrther  fi-om  the  suction  surface.  This  indicates  that  the  code  is  predicting  the 
vortex  slightly  closer  to  the  blade  than  is  observed  experimentally.  The  computed 
pitchwise  flow  angle  distributions  at  80  and  90  percent  pitch  show  a  smaller  leakage  jet 
velocity  and  this,  too,  indicates  a  vortex  located  closer  to  the  blade. 

Figures  2.3(b)  and  2.4(b)  show  the  measured  and  predicted  pitchwise  and 
spanwise  flow  angles,  respectively,  at  an  axial  location  of  80%  chord  from  the  leading 
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edge.  The  spanwise  flow  angle  distributions  show  that  the  vortex  placement  is  still 
nearer  to  the  suction  siirface  for  the  computed  flow  than  in  the  measured  flow.  At  60% 
pitch,  the  measmed  leakage  flow  has  begun  to  wrap  around  back  toward  the  suction 
surface,  while  the  predicted  case  does  not  exhibit  this  trend  until  about  70%  pitch.  At 
80%  pitch,  figure  2.4(b)  shows  the  measured  vortex  to  be  on  its  upward  return  while  the 
predicted  vortex  is  still  penetrating  toward  the  hub.  Figure  2.3(b)  shows  the  same  trends. 
At  70%  pitch,  the  measured  leakage  velocity  is  quite  prevalent,  while  the  predicted  case 
still  shows  signs  of  the  passage  vortex,  with  a  pitchwise  flow  angle  less  than  the  value  far 
from  the  casing.  For  the  region  of  the  passage  closer  to  the  pressure  surface,  the  flow  is 
well  predicted.  This  is  also  the  case  with  the  flow  adjacent  to  the  suction  surface  (90% 
pitch). 

The  flow  angles  at  90%  chord  downstream  of  the  leading  edge  are  presented  in 
part  (c)  of  figures  2.3  and  2.4,  and  the  trends  are  very  much  the  same  as  at  80%  chord. 
The  predicted  vortex  is  still  slightly  closer  to  the  suction  surface  than  the  measured 
vortex.  This  is  most  noticeable  at  60%  pitch  for  the  pitchwise  flow  angle  and  at  70% 
pitch  for  the  spanwise  flow  angle.  Both  measured  and  predicted  angles  begin  to  show  the 
absence  of  the  passage  vortex,  as  it  is  beginning  to  be  pulled  imdemeath  the  leakage 
vortex. 

Overall,  figures  2.3  and  2.4  show  that  the  CFD  code  being  used  is  capable  of 
predicting  the  fine  details  of  the  complex  flow  phenomenon  in  turbomachinery.  The  use 
of  an  isotropic  turbulence  model  in  the  prediction  could  be  accountable  for  the  slight 
discrepancy  in  the  location  of  the  vortex  for  the  predicted  and  measured  flows.  Also, 
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there  is  some  imcertainty  by  the  author  as  to  whether  the  experimental  data  was  taken 
at  pitch  locations  defined  by  the  distance  from  blade  to  blade  or  from  suction  surface  to 
pressure  surface.  If  the  former  case  were  true,  then  this  would  place  the  experimental 
vortex  closer  to  the  suction  surface  and  in  better  agreement  with  the  predicted  case. 
Overall,  the  code  does  a  quality  job  of  predicting  complex  3D  flow  features  and  should 
lend  itself  well  to  the  simulations  presented  in  this  study. 


Spanwise  Distance  Spanwise  Distance  Spanwise  Distance  Spanwise  Distance 
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-  2.5%  Chord  =  Tip  Clearance 

- 1 .0%  Chord  =  Tip  Clearance 

X  2.5%  Chord  =  Tip  Clearance:  Data  Measured  by  Bindon  (1 990) 

+  1 .0%  Chord  =  Tip  Clearance:  Modified  Inlet  Boundary  Conditions  (see  section  4.1) 

(1)  positive  pitch  angle  is  toward  the  pressure  surface. 

(2)  Spanwise  Distance  =  spanwise  distance  from  casing/clearance  height 

(3)  percent  pitch  is  the  percent  of  the  distance  from  pressure 
surface  to  suction  surface. 


Figure  2.3:  Pitch  Angle  near  the  Casing 
Part  (a):  X  =  60%  Axial  Chord 
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-  2.5%  Chord  =  Tip  Clearance 

-  1 .0%  Chord  =  Tip  Clearance 

X  2.5%  Chord  =  Tip  Clearance:  Data  Measured  by  Bindon  (1 990) 

+  1 .0%  Chord  =  Tip  Clearance:  Modified  Inlet  Boundary  Conditions  (see  section  4.1) 

(1)  positive  pitch  angle  is  toward  the  pressure  surface. 

(2)  Spanwise  Distance  =  spanwise  distance  from  casing/clearance  height 

(3)  percent  pitch  is  the  percent  of  the  distance  from  pressure 
surface  to  suction  surface. 


Figure  2.3  (continued):  Pitch  Angle  near  the  Casing 
Part  (a):  X  =  60%  Axial  Chord 
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-  2.5%  Chord  =  Tip  Clearance 

-  1 .0%  Chord  =  Tip  Clearance 

X  2.5%  Chord  =  Tip  Clearance:  Data  Measured  by  Bindon  (1 990) 

+  1 .0%  Chord  =  Tip  Clearance:  Modified  Inlet  Boundary  Conditions  (see  section  4.1) 

(1)  positive  pitch  angle  is  toward  the  pressure  surface. 

(2)  Spanwise  Distance  =  spanwise  distance  from  casing/clearance  height 

(3)  percent  pitch  is  the  percent  of  the  distance  from  pressure 
surface  to  suction  surface. 


Figure  2.3:  Pitch  Angle  near  the  Casing 
Part  (b):  X  =  80%  Axial  Chord 
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-  2.5%  Chord  =  Tip  Clearance 

-  1 .0%  Chord  =  Tip  Clearance 

X  2.5%  Chord  =  Tip  Clearance:  Data  Measured  by  BIndon  (1 990) 

+  1 .0%  Chord  =  Tip  Clearance:  Modified  Inlet  Boundary  Conditions  (see  section  4.1 ) 

(1)  positive  pitch  angle  is  toward  the  pressure  surface. 

(2)  spanwise  Distance  =  spanwise  distance  from  casing/clearance  height 

(3)  percent  pitch  is  the  percent  of  the  distance  from  pressure 
surface  to  suction  surface. 


Figure  2.3  (continued):  Pitch  Angle  near  the  Casing 
Part  (b):  X  =  80%  Axial  Chord 
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-  2.5%  Chord  =  Tip  Clearance 

-  1 .0%  Chord  =  Tip  Clearance 

X  2.5%  Chord  =  Tip  Clearance:  Data  Measured  by  Bindon  (1990) 

+  1 .0%  Chord  =  Tip  Clearance:  Modified  Inlet  Boundary  Conditions  (see  section  4.1 ) 


(1)  positive  pitch  angle  is  toward  the  pressure  surface. 

(2)  Spanwise  Distance  =  spanwise  distance  from  casing/clearance  height 

(3)  percent  pitch  is  the  percent  of  the  distance  from  pressure 
surface  to  suction  surface. 


Figure  2.3  (continued):  Pitch  Angle  near  the  Casing 
Part  (c):  X  =  90%  Axial  Chord 
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-  2.5%  Chord  =  Tip  Clearance 

-  1 .0%  Chord  =  Tip  Clearance 

X  2.5%  Chord  =  Tip  Clearance:  Data  Measured  by  Bindon  (1 990) 

+  1 .0%  Chord  =  Tip  Clearance:  Modified  Inlet  Boundary  Conditions  (see  section  4.1 ) 

(1 )  positive  pitch  angle  is  toward  the  pressure  surface. 

(2)  Spanwise  Distance  =  spanwise  distance  from  casing/clearance  height 

(3)  percent  pitch  is  the  percent  of  the  distance  from  pressure 
surface  to  suction  surface. 


Figure  2.4:  Span  Angle  near  the  Casing 
Part  (a):  X  =  60%  Axial  Chord 
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-  2.5%  Chord  =  Tip  Clearance 

-  1 .0%  Chord  =  Tip  Clearance 

X  2.5%  Chord  =  Tip  Clearance:  Data  Measured  by  Bindon  (1 990) 

+  1 .0%  Chord  =  Tip  Clearance:  Modified  Inlet  Boundary  Conditions  (see  section  4.1 ) 

(1)  positive  pitch  angle  is  toward  the  pressure  surface. 

(2)  Spanwise  Distance  =  spanwise  distance  from  casing/clearance  height 

(3)  percent  pitch  is  the  percent  of  the  distance  from  pressure 
surface  to  suction  surface. 


Figure  2.4  (continued):  Span  Angle  near  the  Casing 
Part  (a):  X  =  60%  Axial  Chord 
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-  2.5%  Chord  =  Tip  Clearance 

- 1 .0%  Chord  =  Tip  Clearance 

X  2.5%  Chord  =  Tip  Clearance:  Data  Measured  by  Bindon  (1 990) 

+  1 .0%  Chord  =  Tip  Clearance:  Modified  Inlet  Boundary  Conditions  (see  section  4.1 ) 

(1 )  positive  pitch  angle  is  toward  the  pressure  surface. 

(2)  Spanwise  Distance  =  spanwise  distance  from  casing/clearance  height 

(3)  percent  pitch  is  the  percent  of  the  distance  from  pressure 
surface  to  suction  surface. 


Figure  2.4:  Span  Angle  near  the  Casing 
Part  (b):  X  =  80%  Axial  Chord 
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-  2.5%  Chord  =  Tip  Clearance 

- 1 .0%  Chord  =  Tip  Clearance 

X  2.5%  Chord  =  Tip  Clearance:  Data  Measured  by  Bindon  (1 990) 

+  1 .0%  Chord  =  Tip  Clearance:  Modified  Inlet  Boundary  Conditions  (see  section  4.1) 

(1 )  positive  pitch  angle  is  towrard  the  pressure  surface. 

(2)  Spanwise  Distance  =  spanwise  distance  from  casing/clearance  height 

(3)  percent  pitch  is  the  percent  of  the  distance  from  pressure 
surface  to  suction  surface. 


Figure  2.4  (continued):  Span  Angle  near  the  Casing 
Part  (b):  X  =  80%  Axial  Chord 
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-  2.5%  Chord  =  Tip  Clearance 

- 1 .0%  Chord  =  Tip  Clearance 

X  2.5%  Chord  =  Tip  Clearance:  Data  Measured  by  Bindon  (1 990) 

+  1 .0%  Chord  =  Tip  Clearance:  Modified  Inlet  Boundary  Conditions  (see  section  4.1 ) 

(1)  positive  pitch  angle  is  toward  the  pressure  surface. 

(2)  Spanwise  Distance  =  spanwise  distance  from  casing/clearance  height 

(3)  percent  pitch  is  the  percent  of  the  distance  from  pressure 
surface  to  suction  surface. 


Figure  2.4  (continued):  Span  Angle  near  the  Casing 
Part  (c):  X  =  90%  Axial  Chord 


CHAPTER  3:  RESULTS 


A  baseline  flat  tip  case  and  the  three  modified  tip  cases  described  in  chapter  1 
were  simulated  using  the  numerical  procedure  outlined  in  chapter  2.  In  this  chapter,  the 
four  computed  flow  fields  are  analyzed  in  detail  and  compared  with  one  another. 

As  a  means  of  qualitatively  determining  the  effect  of  the  modified  blade  tip 
geometries  on  the  leakage  vortex,  particle  traces  through  the  gap  and  leakage  vortex  were 
produced.  Figure  3.1  shows  the  leakage  vortices  for  the  baseline,  leading  edge 
chamfered,  full  chamfered  and  trailing  edge  chamfered  cases.  The  colors  in  Fig.  3.1  are 
added  to  show  similar  regions  of  the  flow  field,  but  should  not  be  considered  to 
correspond  exactly  between  the  various  parts  of  Fig.  3.1 . 

Figure  3.1(a)  shows  the  leakage  vortex  for  the  baseline  case,  which  is  also 
pictured  in  Fig.  1.2.  The  leakage  vortex  develops  along  the  suction  side  comer  of  the 
blade,  beginning  at  about  40%  axial  chord.  All  of  the  fluid  inside  the  leakage  vortex 
passes  through  the  gap  region  with  very  close  proximity  to  the  blade  tip:  roughly  within 
10%  of  the  gap  height.  For  comparison,  the  yellow  fluid  pathlines  in  Fig.  3.1(a)  pass 
through  the  gap  further  removed  from  the  blade  tip,  and  thus,  do  not  roll  into  the  leakage 
vortex.  As  was  discussed  in  the  Background  section,  this  indicates  that  the  rotation  of  the 
leakage  flow  due  to  its  shearing  on  the  blade  tip  is  the  driving  mechanism  for  the 
formation  of  the  leakage  vortex.  The  direction  of  this  rotation  is  opposite  of  the  blade’s 
camber  direction,  and  therefore  negative. 

The  black  fluid  pathlines  in  Fig.  3.1(a)  are  slightly  more  puzzling.  These 
pathlines  also  pass  thorough  the  gap  very  near  to  the  blade  tip,  yet  indicate  an  opposite 
(positive)  rotation  once  inside  the  passage.  The  same  positive  rotation  region  is  shown  in 
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Fig.  1.2  Lab.  iv.  In  the  precursor  study  to  this  one,  Talhnan  and  Lakshminarayana 
[10,11]  discussed  how  the  positive  vorticity  of  the  near-casing  secondary  flow  (Fig.  1.2, 
Lab.  v)  is  carried  into  the  region  below  the  leakage  vortex,  where  it  manifests  itself  as  a 
mildly  rotational  region  of  the  passage  flow  (Fig.  1.2,  Lab.  iv).  To  refer  to  this  region  as 
a  “passage  vortex”  is  incorrect,  however,  since  the  mechanisms  responsible  are  entirely 
different  and  die  rotation  rate  is  considerably  less  than  that  of  a  true  vortex,  such  as  the 
leakage  vortex.  In  comparing  Figs.  1.2  and  3.1(a),  it  appears  that  the  black  fluid 
pathlines  in  Fig.  3.1(a)  are  being  carried  along  with  the  near-casing  secondary  flow  as  it 
is  scraped  below  the  leakage  vortex. 

Figure  3.1(b)  shows  fluid  pathlines  that  describe  the  flow  for  the  leading  edge 
chamfered  case.  The  black  and  blue  fluid  pathlines  illustrate  the  near-casing  secondary 
flow.  As  was  the  design  intention,  this  fluid  is  now  carried  into  the  gap’s  leading  edge 
region.  Comparing  Figs.  3.1(a)  and  3.1(b)  show  that  this  does,  indeed,  act  to  turn  the 
leakage  flow  more  toward  the  camber  direction.  Unfortunately,  the  leakage  flow  has  still 
found  a  way  to  roll-up.  As  has  been  discussed,  the  leakage  flow  in  close  proximity  to  the 
blade  tip  is  responsible  for  causing  the  leakage  vortex.  Figure  3.1(b)  shows  that  for  the 
leading  edge  chamfer  case,  the  red  fluid  pathlines,  while  turned  more  toward  the  camber 
direction,  still  travel  towards  the  suction  side  edge  of  the  blade  underneath  the  obstructing 
near-casing  secondary  flow.  Upon  reaching  the  suction  side  edge,  this  fluid  will  still  roll 
into  a  leakage  vortex. 

For  the  leading  edge  chamfer  case.  Fig.  3.1(b)  shows  that  the  spanAvise  location  of 
the  suction  side  blade  edge  will  vary  with  location.  The  spanwise  location  of  the  suction 
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side  blade  edge  decreases  from  the  leading  edge  to  about  40%  axial  chord,  then  rises 
dramatically  back  to  its  uniform  downstream  value.  This  dramatic  rise  leads  to  an 
enhancement  of  the  leakage  vortex,  as  Fig.  3.1(b)  shows.  The  red  fluid  pathlines  in  Fig. 
3.1(b)  exit  the  gap  region  at  successively  increasing  spanwise  locations.  Therefore,  each 
subsequent  red  fluid  pathline  will  wrap  itself  around  the  leakage  vortex  at  an  increasing 
radius,  leading  to  the  larger  size  of  the  leakage  vortex.  Alternatively,  for  the  baseline 
case  (Fig.  3.1(a)),  the  rolling  up  fluid  pathlines  tend  to  helix  themselves  around  each 
other  at  a  more  constant  radius.  The  larger  leakage  vortex  for  the  leading  edge  chamfer 
case  will  lead  to  an  increase  in  mixing  and  losses  inside  the  passage. 

A  similar  situation  can  be  seen  in  Fig.  3.1(c)  for  the  full  chamfer  case.  Black  and 
blue  fluid  pathlines  indicate  the  obstructing  near-casing  secondary  flow  inside  the  gap, 
while  red  and  yellow  fluid  pathlines  show  the  leakage  flow  passing  underneath.  In  the 
full  chamfer  case,  the  increase  in  spanwise  location  of  the  suction  side  edge  is  spread  out 
more  gradually  across  the  entire  length  of  the  blade.  Therefore,  the  “fluid  stacking” 
effect  seen  in  the  leading  edge  chamfer  case  is  less  severe,  but  still  present.  Once  again, 
the  leakage  vortex  size  and  mixing  is  enhanced  from  that  of  the  baseline  case  in  Fig. 
3.1(a). 

The  flow  through  the  trailing  edge  chamfer  case  is  shown  in  Fig.  3.1(d).  Beyond 
30%  axial  chord,  the  size  of  the  leakage  vortex  is  noticeably  smaller  for  the  trailing  edge 
chamfer  case.  The  reason  for  this  size  reduction  of  the  leakage  vortex  can  be  seen  in  the 
red  fluid  pathlines  of  Fig.  3.1(d).  These  fluid  pathlines  are  the  ones  that  pass  through  the 
gap  near  to  the  blade  tip  and  are  responsible  for  causing  the  leakage  vortex.  Figure  3.1(d) 
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shows  that  from  the  start  of  the  chamfered  region  to  the  blade  trailing  edge,  the  red 
fluid  pathlines  have  been  turned  considerably  towards  the  camber  line’s  direction.  The 
red  fluid  pathlines  in  Fig.  3.1(d)  are  still  sheared  across  the  tip  of  the  blade  and  are,  thus, 
still  rotational.  However,  the  direction  of  this  rotation  is  no  longer  well  aligned  with  the 
direction  of  the  leakage  vortex’s  rotation.  Therefore,  between  40%  and  70%  axial  chord, 
the  initial  leakage  vortex  has  a  more  difficult  time  starting  out.  Further  downstream, 
however,  the  suction  side  edge  of  the  blade  is  increasing  in  its  spanwise  location,  as  it  did 
in  the  leading  edge  and  full  chamfer  cases.  Therefore,  the  leakage  vortex  increases  in 
size  through  the  final  30%  of  the  passage. 

Because  only  the  trailing  edge  chamfer  case  shows  a  reduction  in  the  leakage 
vortex,  the  remainder  of  this  report  will  focus  on  comparison  between  the  baseline  and 
trailing  edge  chamfer  cases.  A  comparison  of  the  flow  in  the  gap  region  for  these  two 
cases  is  given  in  Fig.  3.2.  Once  again,  fluid  pathlines  are  used  to  illustrate  the  flow 
qualitatively.  Figure  3.2  clearly  shows  the  turning  of  the  leakage  flow  towards  the 
camber  direction  that  results  from  the  blade  tip  modification.  The  red  fluid  pathlines  in 
the  trailing  edge  chamfer  case  (Fig.  3.2(b))  are  at  an  angle  of  approximately  45  degrees 
from  the  camber  direction,  while  in  the  baseline  case  (Fig.  3.2(a)),  they  are  almost  normal 
to  the  camber  direction.  In  the  upstream  half  of  the  gap,  however,  the  two  leakage  flows 
are  very  similar.  The  turning  of  the  leakage  flow  near  the  blade  tip  in  the  trailing  edge 
chamfer  case  is  of  great  significance  in  attempts  to  reduce  the  leakage  vortex,  since  a 
leakage  flow  aligned  with  the  camber  direction  cannot  cause  a  leakage  vortex.  Also,  the 
red  and  yellow  fluid  pathlines  in  Fig.  3.2(a)  show  that  the  leakage  flow  in  the  trailing 


51 


edge  chamfer  case  is  highly  three-dimensional.  The  yellow  fluid  pathlines  in  Fig. 
3.2(b)  are  only  slightly  further  removed  from  the  blade  tip  than  the  red  fluid  pathlines. 

In  both  parts  of  Fig.  3.2,  the  black  fluid  pathlines  illustrate  the  near-casing 
secondary  flow.  These  fluid  pathlines  are  released  from  the  same  starting  point  for  both 
the  baseline  case  and  the  trailing  edge  chamfer  case.  In  the  baseline  case,  some  near¬ 
casing  secondary  flow  enters  the  gap  region  near  the  leading  edge  of  the  blade,  but  this 
fluid  is  forced  back  out  of  the  gap  by  the  leakage  flow  only  slightly  further  downstream 
from  where  it  entered.  In  the  trailing  edge  chamfer  case.  Fig.  3.2(b)  shows  that  the  near¬ 
casing  secondary  flow  entering  the  gap  near  the  leading  edge  is  able  to  continue  inside 
the  gap  longer,  with  a  direction  more  similar  to  the  blade’s  camber.  This  is  due  to  the 
increased  flow  area  available  in  the  downstream  half  of  the  gap.  In  addition,  near-casing 
secondary  flow  is  able  to  enter  the  gap  region  further  downstream  in  the  trailing  edge 
chamfer  case  than  in  the  baseline  case. 

The  increased  quantity  of  near-casing  secondary  flow  and  increased  longevity  of 
said  flow  inside  the  gap  region  are  responsible  for  the  turning  of  the  leakage  flow  toward 
the  camber  direction.  The  near-casing  secondary  flow  that  enters  the  gap  acts  as  an 
obstruction  to  the  leakage  flow.  The  leakage  flow  then  has  two  ways  that  it  can 
circumnavigate  the  obstruction.  First,  the  leakage  flow  can  force  its  way  underneath  the 
obstruction,  such  as  what  is  seen  in  the  leading  edge  region  of  Fig.  3.2(b).  Second,  the 
leakage  flow  can  turn  toward  the  camber  direction  and  pass  around  the  obstruction 
further  downstream.  Figure  3.2(b)  indicates  that  the  leakage  flow  utilizes  both  of  these 
options  in  the  trailing  edge  chamfer  case.  The  turning  of  the  yellow  and  red  leakage  flow 
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pathlines  indicates  the  latter  option,  while  the  crossing  of  the  black  and  yellow  fluid 
pathlines  indicates  the  former. 

Figure  3.3  shows  both  the  velocity  and  pressure  field  in  the  blade-to-blade  plane 
at  a  spanwise  location  of  100.2%  of  the  blade  hei^t.  This  location  is  20%  of  the 
distance  fi-om  the  top  of  the  blade  to  the  outer  casing.  Once  again,  a  direct  comparison 
between  the  baseline  and  trailing  edge  chamfer  cases  is  made.  Pressure  is  displayed  as 
pressure  coefficient  (Cp)  contours,  where  the  pressure  coefficient  is  defined  as  the 
difference  between  the  local  and  inlet  static  pressure  values,  non-dimensionalized  by  the 
inlet  dynamic  pressure.  A  line  has  been  added  to  Fig.  3.3(b)  to  indicate  the  location 
where  the  trailing  edge  chamfering  begins.  The  same  line  has  been  added  to  Fig.  3.3(a) 
as  a  reference. 

For  the  baseline  case.  Fig.  3.3(a)  shows  that  the  velocity  and  pressure  fields  inside 
the  gap  are  well  aligned  with  each  other.  Velocity  vectors  point  in  a  direction  of 
decreasing  pressure,  and  that  direction  is  fairly  uniform  throughout  the  gap  region.  The 
leakage  flow  exits  the  gap  in  a  direction  fairly  normal  to  the  suction  surface  of  the  blade. 
Therefore,  the  vorticity  of  the  leakage  flow  near  the  blade  tip  will  be  similar  to  the 
direction  of  the  leakage  vortex’s  rotation.  For  the  trailing  edge  chamfer  case.  Fig.  3.3(b) 
shows  velocity  vectors  from  the  leakage  flow  and  near-casing  secondary  flow  turmng  on 
each  other  at  the  suction  side  of  the  gap.  The  leakage  flow  velocity  vectors  turn  toward 
the  camber  direction,  as  was  seen  in  Figs.  3.1(d)  and  3.2(b).  Caution  should  be  taken 
when  following  the  velocity  vectors  through  the  gap  in  Fig.  3.3(b),  however,  since  Figs. 
3.1(d)  and  3.2(b)  showed  an  important  spanwise  component  to  the  flow. 
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The  original  reasoning  behind  the  trailing  edge  chamfer  design  was  that  the 
increased  area  made  available  to  the  leakage  flow  would  act  to  diffuse  the  flow  and 
reduce  the  leakage  velocity.  The  fluid  pathlines  in  Figs.  3.1(d)  and  3.2(b),  however, 
indicate  that  the  flow  inside  the  chamfered  gap  is  considerably  more  complex  than  that 
through  a  simple  ID  diffuser.  A  pressure  rise  in  the  direction  of  the  leakage  flow  is  seen 
in  the  trailing  edge  region  of  Fig.  3.3(b),  which  would  indicate  that  the  leakage  flow  is 
being  decelerated.  Comparison  with  Fig.  3.3(a)  indicates  that  the  leakage  flow  is  indeed 
loosing  momentum  in  the  direction  corresponding  to  the  leakage  flow  in  the  baseline 
case.  The  leakage  flow  does  not  decelerate,  however,  and  instead  turns  toward  the 
trailing  edge,  where  lower  pressure  values  exist.  Therefore,  the  turning  of  the  leakage 
flow  toward  the  camber  direction  in  the  trailing  edge  chamfer  case  is  due  to  two 
mechanisms:  1)  the  obstruction  of  the  leakage  flow  by  the  near-casing  secondary  flow 
and  2)  the  pressure  rise  across  the  chamfered  gap  due  to  its  diffusing  effect.  With  these 
two  mechanisms,  the  pressure  rise  in  the  latter  mechanism  may  be  entirely  due  to  the 
former  mechanism  and  not  a  diffuser  effect  at  all.  Throughout  Fig.  3.3(b),  the  velocity 
and  pressure  fields  appear  to  be  well  coupled  together.  That  is,  velocity  vectors  tend  in 
direction  toward  lower  Cp  values. 

The  difference  in  the  direction  of  the  red  and  yellow  fluid  pathlines  of  Fig.  3.2(b) 
can  also  be  explained  by  considering  the  pressure  field  inside  the  gap.  Figures  similar  to 
Fig.  3.3(b)  (not  shown)  indicate  that  the  pressure  distribution  in  the  blade-to-blade  plane 
does  not  change  significantly  inside  the  gap  with  spanwise  location.  Therefore,  the 
increased  turning  of  the  red  fluid  pathlines  in  Fig.  3.2(b)  is  similar  to  a  passage  vortex 
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effect.  The  flow  in  the  gap  establishes  a  pressure  gradient  to  support  the  turning  of  the 
leakage  flow  toward  the  camber  direction.  At  locations  closer  to  the  blade  tip,  the  fluid  is 
slower  moving  and,  thus,  requires  less  centripetal  acceleration  than  the  pressure  gradient 
is  providing.  Therefore,  the  fluid  is  overturned,  just  as  a  boundary  layer  in  a  curved  duct 
becomes  overturned.  This  overturning  of  the  leakage  flow  near  the  blade  tip  is  beneficial, 
since  this  fluid  is  exactly  the  fluid  that  forms  the  leakage  vortex. 

Figure  3.4  is  similar  to  Fig.  3.3,  only  at  a  spanwise  location  of  98%  blade  span. 
The  darkest  blue  contour  indicates  the  core  of  the  leakage  vortex  in  both  parts.  This  core 
is  noticeably  smaller  in  the  trailing  edge  chamfer  case.  Furthermore,  the  beginning  of  the 
core  has  moved  toward  the  trailing  edge  of  the  blade  with  the  trailing  edge  chamfer  case. 
The  inception  of  the  leakage  vortex  has  been  delayed  with  the  modified  blade  tip. 
Immediately  adjacent  to  the  leakage  vortex  core,  the  trailing  edge  chamfer  case  shows 
reduced  pressure  values,  which  closer  resemble  the  pressure  contours  at  the  blade’s  mid¬ 
span  location  (not  shown).  This  indicates  improved  turbine  performance  near  the  tip  for 
the  trailing  edge  chamfer  case.  Upstream  and  throughout  the  remainder  of  the  passage, 
the  flow  fields  are  similar  for  the  two  cases. 

Figures  3.5  through  3.8  and  Fig.  3.10  show  secondary  vorticity  contours  in  the 
tangential-spanwise  plane  near  the  leakage  vortex  at  various  axial  locations.  Vorticity 
shown  is  in  the  2D-streamwise  direction,  and  is  non-dimensionalized  using  the  blade 
span  and  inlet  velocity.  Figure  3.5  is  situated  at  60%  axial  chord.  Both  parts  of  Fig.  3.5 
indicate  that  vorticity  is  being  convected  out  of  the  gap  region,  with  negative  vorticity 
coming  fi'om  near  the  tip  and  positive  vorticity  coming  from  near  the  outer  casing.  The 
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negative  vorticity  region  inside  the  gap  is  noticeably  larger  in  the  trailing  edge  chamfer 
case.  This  is  related  to  the  three  dimensionality  of  the  leakage  flow  in  the  downstream 
half  of  the  trailing  edge  chamfer  case,  illustrated  by  the  yellow  and  red  fluid  pathlines  in 
Fig.  3.2(b).  These  pathlines  show  that  tangential  velocity  decreases  with  proximity  to  the 
blade  tip,  resulting  in  the  increased  levels  of  negative  vorticity  inside  the  gap  for  the 
trailing  edge  chamfer  case. 

In  Fig.  3.5,  the  region  of  negative  vorticity  is  associated  with  the  flow  in  and 
aroimd  the  leakage  vortex  in  both  cases.  Above  and  to  the  right  of  this  region  is  a 
complicated  region  of  positive  vorticity,  associated  with  the  near-casing  secondary  flow. 
Above  the  leakage  vortex,  the  region  of  greatest  positive  vorticity  is  associated  with  the 
turning  of  the  near-casing  secondary  flow  away  fi'om  the  gap.  This  turning  was 
illustrated  by  the  black  fluid  pathlines  in  Fig.  3.2.  With  the  trailing  edge  chamfer  case, 
the  turning  of  the  near-casing  secondary  flow  takes  place  closer  to  the  gap,  where  it 
further  inhibits  the  leakage  flow.  The  region  of  slightly  less  positive  vorticity  to  the  right 
of  the  leakage  vortex  is  composed  of  near-casing  secondary  flow  and  leakage  flow  from 
near  the  casing  that  have  exited  the  suction  side  of  the  gap  further  upstream.  Because  of 
the  casing’s  relative  motion,  both  of  these  types  of  fluid  will  be  positively  rotational,  and 
that  rotation  later  appears  inside  the  passage  as  the  region  aforementioned. 

Figures  3.6  through  3.8  are  identical  to  Fig.  3.5,  only  at  axial  locations  of  70%, 
80%,  and  90%  axial  chord,  respectively.  All  of  the  trends  discussed  in  Fig.  3.5  are  also 
present  in  Figs.  3.6  through  3.8.  Together,  the  four  figures  illustrate  the  development  of 
the  secondary  flow  field  in  the  downstream  half  of  the  passage.  The  four  figures  show 
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that  the  growth  of  the  region  of  negative  vorticity  associated  with  the  leakage  vortex  is 
less  severe  in  the  trailing  edge  chamfer  case.  By  90%  axial  chord,  the  size  of  this  region 
in  the  trailing  edge  chamfer  case  is  approximately  half  of  its  value  in  the  baseline  case. 
This  agrees  with  the  illustrations  of  the  leakage  vortex  that  were  shown  in  Figs.  3.1  and 
3.2.  Additionally,  the  region  of  positive  vorticity  below  and  to  the  right  of  the  leakage 
vortex  shows  reduced  values  for  the  trailing  edge  case  in  all  of  Figs.  3.5  through  3.8, 
indicating  improved  flow  performance  inside  the  passage. 

Secondary  velocity  vectors  are  shown  along  with  the  vorticity  contours  in  Fig. 
3.8.  At  first  glance,  the  velocity  vectors  inside  the  gap  appear  reduced  in  velocity 
magnitude  for  the  trailing  edge  case.  However,  care  must  be  taken  in  analyzing  2D 
velocity  vectors,  as  a  third  component  is  neglected.  Instead,  close  analysis  of  Fig.  3.2(b) 
indicates  that  the  flow  inside  the  gap  is  not  decelerated,  but  rather  turned  toward  the 
camber  direction.  This  turning  was  seen  in  Figs.  3.1  and  3.2.  Leakage  vortex  roll-up  is 
much  more  clearly  observable  in  the  velocity  vectors  for  the  baseline  case.  The  baseline 
case  in  Fig.  3.8  shows  a  stronger  flow  toward  the  outer  casing  on  the  left  side  of  the 
leakage  vortex.  The  shearing  of  this  secondary  flow  by  the  blade’s  suction  surface  results 
in  the  region  of  positive  vorticity  to  the  left  of  the  leakage  vortex.  The  shearing  results  in 
increased  profile  losses  and  is,  thus,  undesirable.  This  secondary  flow  is  considerably 
reduced  for  the  trailing  edge  chamfer  case.  Therefore,  the  trailing  edge  chamfer  case 
once  again  shows  improved  performance  over  the  baseline  case.  Additionally,  Fig.  3.9 
shows  the  actual  secondary  vorticity  values  inside  the  leakage  vortex  of  Fig.  3.8.  Figure 
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3.9  shows  quantitatively  that  the  secondary  vorticity  magnitude  at  the  core  of  the 
vortex  has  reduced  for  the  trailing  edge  chamfer  case. 

Figure  3.10  shows  one  final  plot  of  the  vorticity  field:  this  time  downstream  of  the 
trailing  edge  at  110%  axial  chord.  This  figure  clearly  indicates  that  the  effect  of  the 
trailing  edge  chamfer  is  to  reduce  the  strength  of  the  leakage  vortex.  Downstream  of  the 
passage,  the  contours  of  negative  vorticity  associated  with  the  leakage  vortex  are  reduced 
in  both  size  and  magnitude  with  the  trailing  edge  chamfer  case.  This  will  result  in 
reduced  downstream  losses  associated  with  the  decay  of  the  leakage  vortex,  as  well  as 
reduced  relative  unsteadiness  in  subsequent  blade  rows.  Additionally,  the  reduction  of 
the  leakage  vortex’s  vorticity  seen  in  Fig.  3.10  could  be  partially  due  to  increased 
turbulence  production  inside  the  gap,  although  analysis  of  the  turbulence  field  could  not 
confirm  or  deny  this  possibility.  Figure  3.11  shows  the  secondary  vorticity  values  inside 
the  center  of  the  leakage  vortices  of  Fig.  3.10.  Figure  3.1 1  once  again  shows  that  the  size 
and  strength  of  the  leakage  vortex  is  reduced  with  trailing  edge  chamfering. 

Figures  3.12  and  3.13  show  contours  of  total  pressure  in  the  tangential-spanwise 
plane.  In  these  two  figures,  total  pressure  is  normalized  by  the  inlet  passage-averaged 
total  pressure  plus  an  additional  work  component.  The  work  component  is  the  additional 
total  pressure  made  available  to  the  flow  via  the  work  done  on  the  flow  by  the  relatively 
moving  outer  casing.  This  work  term  increases  with  axial  location,  but  does  not  exceed 
more  than  3%  of  the  inlet  total  pressure. 

Figure  3.12  shows  total  pressure  contours  at  90%  axial  chord.  Initially,  Fig.  3.12 
appears  to  contradict  the  findings  presented  earlier.  The  region  of  lowest  total  pressure  is 
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situated  around  the  leakage  vortex.  This  region  is  larger  in  the  trailing  edge  chamfer 
case,  despite  the  case’s  smaller  leakage  vortex.  To  explain  this  feature,  the  total  pressure 
throughout  the  flow  needs  to  be  analyzed  in  greater  detail.  Figure  3.12  shows  higher 
losses  inside  the  gap  for  the  trailing  edge  chamfer  case.  A  closer  analysis  of  the  total 
pressure  field  shows  that  this  trend  is  consistent  throughout  the  chamfered  region  of  the 
gap.  Figures  showing  this  trend  are  not  presented,  but  evidence  of  the  trend  is  present  in 
Fig.  3.2.  hi  Fig.  3.2(b),  the  red  and  yellow  fluid  pathlines  exhibit  shearing  across  one 
another,  which  will  lead  to  viscous  dissipation  of  energy  and  losses.  The  largest  total  loss 
contours  in  Fig.  3.12  are  therefore  not  associated  entirely  with  the  leakage  vortex,  but 
with  the  leakage  vortex  and  the  high  loss  fluid  convected  out  of  the  gap  along  pathlines 
such  as  the  yellow  ones  in  Fig.  3.1(d).  Figure  3.12  shows  that  the  losses  away  from  the 
leakage  vortex,  associated  with  the  blade’s  boundary  layer,  are  comparable  for  the  two 
cases.  Therefore,  at  90%  axial  chord,  the  trailing  edge  chamfer  case  does  not  show  an 
improvement  in  losses,  despite  the  smaller  leakage  vortex.  Additionally,  Fig.  3.13  shows 
total  pressure  values  inside  the  leakage  vortices  shown  in  Fig.  3.12.  For  the  trailing  edge 
chamfer  case.  Fig.  3.13  shows  decreased  losses  associated  with  the  core  of  the  leakage 
vortex  and  increased  losses  in  the  region  adjacent  to  the  leakage  vortex. 

Figure  3.14  shows  total  pressure  contours  at  1 10%  axial  chord.  In  both  cases,  the 
contour  of  lowest  total  pressure  is  associated  directly  with  the  decaying  leakage  vortex,  as 
the  leakage  flow  is  no  longer  present  downstream  of  the  blade.  The  losses  associated 
with  the  decaying  leakage  vortex  are  smaller  in  the  trailing  edge  chamfer  case,  as  the  blue 
contours  show.  While  the  region  of  reduced  total  pressure  in  Fig.  3.14  occupies  a 
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similarly  sized  area  for  both  cases,  total  pressure  values  appear  slightly  higher  in  the 
trailing  edge  chamfer  case.  This  agrees  with  the  reduced  overall  secondary  flow  activity 
in  the  trailing  edge  chamfer  case,  which  was  observed  in  previous  figures.  In  the  trailing 
edge  chamfer  case,  increased  losses  have  taken  place  inside  the  gap  further  upstream,  and 
the  resulting  reduced  total  pressure  values  are  more  evenly  spread  out  at  110%  axial 
chord,  hi  the  baseline  case,  losses  are  being  produced  at  110%  axial  chord  to  a  greater 
extent,  due  to  the  more  prevalent  leakage  vortex  and  surrounding  secondary  flows. 
Figure  3.15  shows  the  total  pressure  values  inside  the  leakage  vortices  in  Fig.  3.14.  For 
the  trailing  edge  chamfer  case,  losses  are  again  reduced  inside  the  leakage  vortex  core 
and  increased  in  the  region  surrounding  the  leakage  vortex. 

Figure  3.16  shows  the  passage-averaged  loss  coefficient  values  versus  axial 
location  for  the  four  cases  shown  in  Fig.  3.1.  The  passage-averaged  loss  coefficient  is 
defined  as  the  difference  between  the  inlet  and  local  total  pressure  values,  mass  averaged 
over  the  tangential-spanwise  plane.  This  total  pressure  difference  is  then  non- 
dimensionalized  by  the  inlet  dynamic  pressure.  Please  also  note  that  the  value  of  the  inlet 
total  pressure  is  increased  slightly  with  axial  location,  in  order  to  account  for  the  work 
done  on  the  fluid  by  the  moving  outer  casing  in  the  relative  frame  of  reference. 

The  leading  edge  chamfer  and  full  chamfer  cases  exhibit  higher  losses  throughout 
the  blade  row.  This  is  to  be  expected,  as  both  cases  showed  a  larger  leakage  vortex  than 
the  baseline  case.  For  the  baseline  and  trailing  edge  chamfer  cases,  the  loss  coefficient 
values  inside  the  blade  row  are  of  negligible  difference.  A  difference  in  loss  coefficient 
of  0.01  corresponds  to  a  difference  in  passage  averaged  total  pressure  of  about  1  pascal. 
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which  is  less  than  the  numerical  accuracy  of  the  solution  and  passage-averaging 
integration.  Therefore,  it  appears  that  with  the  trailing  edge  chamfer  modification,  the 
increase  in  losses  inside  the  gap  region  offset  the  decrease  in  losses  due  to  decreased 
secondary  flow  activity  inside  the  passage. 

Downstream  of  the  blade,  the  trailing  edge  chamfer  case  exhibits  lower  losses 
than  the  baseline  case.  The  decrease  in  losses  is  due  to  the  smaller  size  of  the  leakage 
vortex.  Losses  associated  with  the  tip  leakage  vortex  are  due  to  the  mixing  of  the  leakage 
vortex  and  viscous  dissipation.  Downstream  of  the  blade  row,  the  leakage  vortex  is  no 
longer  being  supplied  with  energy  from  the  leakage  flow  and  will  eventually  decay  away. 
The  smaller  vortex  in  the  trailing  edge  chamfer  case  decays  more  rapidly  and  with  less 
integrated  loss  of  energy  than  the  larger  leakage  vortex  in  the  baseline  case. 

The  full  chamfer  case  also  exhibits  reduced  passage-averaged  losses  downstream 
of  the  blade  row,  although  the  difference  in  loss  coefficient  values  between  the  full 
chamfer  and  trailing  edge  chamfer  cases  is  once  again  smaller  than  the  numerical 
accuracy  of  the  calculating  procedure.  The  full  chamfer  case  should  show  increased 
losses  downstream  of  the  blade  row  due  to  the  increased  size  of  the  leakage  vortex. 
Therefore,  the  reduction  in  downstream  losses  for  the  full  chamfer  case  is  likely  due  to  a 
reduction  in  the  losses  away  from  the  leakage  vortex.  A  comparison  of  the  black  fluid 
pathlines  in  Figs.  3.1(a)  and  3.1(c)  supports  this  conclusion,  as  there  is  less  fluid  rotation 
exhibited  in  these  pathlines  for  the  full  chamfer  case.  Still,  passage-averaged  losses  are 
not  the  only  criteria  for  judging  the  performance  of  the  blade.  A  smaller  leakage  vortex 
will  result  in  less  unsteadiness  within  a  subsequent  downstream  stator.  This  will  lead  to  a 
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reduction  in  unsteady  blade  loading  and  noise  downstream.  Therefore,  the  trailing 
edge  chamfer  case  is  still  advantageous  over  the  full  chamfer  case. 

One  further  criterion  for  judging  turbine  blade  performance  is  the  loading  on  the 
blade.  By  reducing  the  turning  of  the  blade  to  zero  at  the  tip,  the  leakage  vortex  is 
eliminated,  but  so  is  the  work  being  done  on  the  blade  by  the  flow.  Therefore,  caution 
has  to  be  taken  with  techniques  that  reduce  the  leakage  vortex. 

Figure  3.17  shows  contours  of  static  pressure  coefficient  on  the  blade  surface  for 
the  two  cases.  Static  pressure  coefficient  (Cp)  was  defined  earlier  in  the  text,  in  reference 
to  Figs.  3.3  and  3.4.  While  the  pressure  surface  of  the  blade  is  not  shown  in  Fig.  3.17,  the 
pressure  distribution  on  that  surface  is  almost  identical  for  the  two  cases.  The  Cp 
distribution  on  the  blade  tip  is  very  similar  to  the  Cp  distribution  inside  the  gap  in  Fig. 
3.3,  indicating  that  pressure  does  not  change  dramatically  in  the  spanwise  direction 
within  the  gap.  The  contour  of  lowest  Cp  values  situated  along  the  suction  side/blade  tip 
comer  is  associated  with  the  leakage  vortex  in  both  cases.  This  contour  is  smaller  in  size 
and  delayed  in  its  origin,  agreeing  with  the  illustrations  of  the  two  leakage  vortices  in  Fig. 
3.1.  Away  fi:om  the  leakage  vortex,  the  Cp  distribution  on  the  suction  surface  is  similar 
for  the  two  cases,  indicating  similar  blade  loading. 

Figure  3.18  shows  the  Cp  distribution  on  the  blade  at  the  tip  for  the  baseline  and 
trailing  edge  chamfer  case.  Additionally,  the  mid-span  Cp  distribution  is  also  shown. 
Before  30%  axial  chord,  the  two  blade  tip  Cp  distributions  are  very  similar.  This  is 
expected,  since  chamfering  does  not  begin  until  further  downstream.  On  the  suction  side 
of  the  blade,  the  Cp  distribution  for  the  trailing  edge  chamfer  case  closer  resembles  the 


62 


distribution  at  the  mid-span.  The  sudden  jump  in  the  trailing  edge  chamfer  case’s  Cp 
distribution  on  the  suction  surface  at  45%  axial  chord  is  associated  with  the  start  of  the 
blade  chamfer.  The  miriimum  Cp  value  on  the  suction  surface  blade  tip  is  of  similar 
value  for  the  two  cases.  On  the  pressure  surface,  Cp  values  are  reduced  from  the  mid¬ 
span  values  for  both  cases,  indicating  the  unloading  of  the  blade  tip  from  the  entrainment 
of  the  leakage  flow  into  the  gap.  With  the  trailing  edge  chamfering,  this  unloading 
increases  slightly  (reduced  Cp)  through  the  upstream  50%  axial  chord  and  decreases 
slightly  (increased  Cp)  through  the  downstream  50%  axial  chord.  In  general,  the  pressure 
drop  across  the  blade  tip  is  smaller  in  the  trailing  edge  chamfer  case. 

The  integrated  forces  acting  on  the  blade  are  shown  in  Table  3.1  for  the  four  cases 
being  considered.  Force  values  are  reported  as  the  ratio  of  a  particular  case’s  value  to  the 
baseline  case’s  value,  as  well  as  a  percent  change  from  the  baseline  case.  The  tangential 
force  component  provides  torque  to  the  rotor  and,  thus,  power  output  from  the  turbine. 
The  power  supplied  to  the  turbine  blade  is  simply  the  blade  loading  in  the  tangential 
direction  times  the  blade’s  velocity  (31.0  m/s  in  this  case).  This  means  that  the 
percentage  change  in  blade  loading  is  equivalent  to  the  percentage  change  in  power  to  the 
blade.  The  tangential  blade  loading  for  the  trailing  edge  chamfer  case  is  reduced  from  the 
baseline  case  by  1.03%.  While  this  number  is  quite  small,  the  additive  effect  over 
multiple  blades  in  a  rotor  could  be  non-negligible.  The  axial  component  of  the  blade 
loading  is  primarily  due  to  the  pressure  drop  through  the  turbine  and  does  not  contribute 
to  the  torque  applied  to  the  turbine’s  axis.  For  the  trailing  edge  chamfer  case,  this 
component  is  reduced  by  1.71%  from  the  baseline  case’s  value.  Depending  on  the 


design,  a  small  cost  in  blade  loading  may  be  an  acceptable  trade-off  for  improved 
noise,  unsteadiness,  durability,  and  downstream  performance  that  comes  with  a  reduced 
leakage  vortex. 

Table  3.1:  Blade  Loading  Force  Values 

Total  Force  /  I  Baseline  Leading  Edge  I  Full  Chamfer  Trailing  Edge 
Baseline  Total  Force  Case  Chamfer  Case  Case  Chamfer  Case 
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Figure  3,5 ;  Non-Dimensional  Streamwise  Vorticity  Contours  at  60%  Axial  Chord 
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Figure  3.6;  Non-Dimensional  Streamwise  Vorticity  Contours  at  70%  Axial  Chord 
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Figure  3.7:  Non-Dimensional  Streamwise  Vorticity  Contours  at  80%  Axial  Chord 
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Figure  3,8:  Non-Dimensional  Streamwise  Vorticity  Contours  at  90%  Axial  Chord 
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Figure  3.9:  Streamwise  Vorticity  Inside  the  Leakage  Vortex  at  90%  Axial  Chord 
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110%  Axial  Chord 
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Figure  3.10:  Non-Dimensional  Streamwise  Vorticity  Contours  at  110%  Axial  Chord 
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Figure  3.11;  Streamwise  Vorticity  Inside  the  Leakage  Vortex  at  110%  Axial  Chord 
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Figure  3.12:  Total  Pressure  /  Inlet  Total  Pressure  at  90%  Axial  Chord 


Total  Pressure  /  Inlet  Total  Pressure  Total  Pressure  /  Inlet  Total  Pressure 


Percent  Blade  Span  Percent  Blade  Span 


110%  Axial  Chord 


Baseline  Case 


0  0.25  0.5  0.75  1 

(Y  -  Yperiodic)  /  Blade  Spacing 


1 1 0%  Axial  Chord  TE  Chamfer  Case 


0  0.25  0.5  0.75  1 

(Y  -  Yperiodic)  /  Blade  Spacing 


Figure  3.14:  Total  Pressure  /  Inlet  Total  Pressure  at  110%  Axial  Chord 
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Figure  3.15:  Total  Pressure  Inside  the  Leakage  Vortex  at  110%  Axial  Chord 
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Figure  3.16;  Passage-Average  Loss  CoefRcient  vs.  Axial  Location 


Figure  3.17:  Static  Pressure  Coefficient  (Cp)  on  the  Blade 


Cp  =  (Pstatic  -  Pinlet)  /  inlet  Dynamic  Pressure 
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Black:  Mid-Span 
Red:  Baseline  Case 


Figure  3.18:  Static  Pressure  Coefficient  (Cp)  at  the  Blade  Tip  and  Midspan 


CHAPTER  4:  SUMMARY  AND  CONCLUSIONS 


Three  different  tip  geometry  modifications  were  applied  to  an  axial  flow  turbine 
case.  These  modified  blade  tips  were  designed  based  on  current  understanding  of  the 
flow  physics  in  axial  flow  turbines,  with  the  intention  of  reducing  the  leakage  flow  and 
vortex.  The  three  designs  involved  adding  a  chamfer  to  the  leading  edge  region  of  the 
tip,  the  trailing  edge  region  of  the  tip,  and  the  entire  tip  region.  The  angle  of  the  chamfer 
was  7.7  degrees  from  horizontal,  and  the  gap  height  was  kept  uniform  along  the  pressure 
side  of  the  gap  (1.0%  blade  span). 

Along  with  a  baseline,  flat  tipped  case,  the  flow  fields  for  the  three  modified  blade 
tip  designs  were  simulated  using  a  pressure-correction  type,  3D  Navier-Stokes  CFD 
solver  with  k-e  turbulence  modeling.  The  solution  flow  fields  were  then  analyzed  in 
detail. 

Chamfering  the  leading  edge  of  the  blade  tip  led  to  a  larger  leakage  vortex  and 
increased  losses.  Chamfering  the  entire  length  of  the  blade  tip  also  led  to  a  larger  leakage 
vortex,  but  showed  potential  for  turning  the  leakage  flow  more  toward  the  camber 
direction.  Chamfering  only  the  trailing  edge  region  of  the  blade  tip  caused  the  leakage 
flow  very  near  to  the  blade  tip  to  turn  toward  the  camber  direction  significantly.  The 
turning  of  the  leakage  flow  was  found  to  be  primarily  due  to  the  obstruction  of  the  near¬ 
casing  secondary  flow  inside  the  gap  near  the  suction  surface.  The  turned  leakage  flow 
was  found  to  exit  the  gap  at  an  angle  less  normal  to  the  direction  of  the  leakage  vortex’s 
rotation.  The  vorticity  of  the  leakage  flow  near  the  blade  tip  was,  therefore,  less  aligned 
with  the  rotation  of  the  leakage  vortex.  The  resulting  leakage  vortex  for  the  trailing  edge 
chamfer  case  was  reduced  in  size  and  of  different  shape  than  the  leakage  vortex  in  the 
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baseline  case.  The  smaller  leakage  vortex  in  the  trailing  edge  chamfer  case  showed 
reduced  losses  downstream  of  the  passage.  The  reduced  size  of  the  vortex  would  also 
lead  to  reduced  imsteadiness  and  three-dimensionality  downstream  and,  thus,  improved 
performance  in  subsequent  blade  rows. 

Chamfering  the  trailing  edge  of  the  blade  tip  has  been  shown  to  reduce  the  size  of 
the  leakage  vortex  in  axial  flow  turbines.  Future  blade  designs  should  take  advantage  of 
the  leakage  flow  turning  that  is  present  in  these  cases.  Future  studies  in  this  area  should 
investigate  the  performance  of  trailing  edge  chamfered  blade  tips  experimentally  to  verify 
their  improved  performance.  In  addition,  the  chamfer  angle  of  7.7  degrees  was  chosen 
somewhat  arbitrarily,  and  additional  studies  should  further  investigate  how  the  choice  of 
chamfer  angle  affects  performance.  Finally,  the  addition  of  a  small  fence  along  the 
suction  side  of  the  blade  tip’s  chamfered  region  may  allow  for  the  turning  of  the  leakage 
flow  while  prohibiting  the  vortex-causing  near-tip  leakage  flow  from  exiting  the  gap. 
The  latter  two  recommendations  for  future  research  would  lend  themselves  well  to  CFD 


analysis. 
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